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Abstract 

Contraction, expansion and despinning have been common in the past evolution of Solar 
System bodies. These processes deform the lithosphere until it breaks along faults. Their 
characteristic tectonic patterns have thus been sought for on all planets and large satellites 
with an ancient surface. While the search for despinning tectonics has not been conclu- 
sive, there is good observational evidence on several bodies for the global faulting pattern 
associated with contraction or expansion, though the pattern is seldom isotropic as pre- 
dicted. The cause of the non-random orientation of the faults has been attributed either 
to regional stresses or to the combined action of contraction/expansion with another de- 
formation (despinning, tidal deformation, reorientation). Another cause of the mismatch 
may be the neglect of the lithospheric thinning at the equator or at the poles due either 
to latitudinal variation in solar insolation or to localized tidal dissipation. Using thin elas- 
tic shells with variable thickness, I show that the equatorial thinning of the lithosphere 
transforms the homogeneous and isotropic fault pattern caused by contraction/expansion 
into a pattern of faults striking east-west, preferably formed in the equatorial region. By 
contrast, lithospheric thickness variations only weakly affect the despinning faulting pattern 
consisting of equatorial strike-slip faults and polar normal faults. If contraction is added to 
despinning, the despinning pattern first shifts to thrust faults striking north-south and then 
to thrust faults striking east-west. If the lithosphere is thinner at the poles, the tectonic 
pattern caused by contraction/expansion consists of faults striking north/south. I start by 
predicting the main characteristics of the stress pattern with symmetry arguments. I fur- 
ther prove that the solutions for contraction and despinning are dual if the inverse elastic 
thickness is limited to harmonic degree two, making it easy to determine fault orientation 
for combined contraction and despinning. I give two methods for solving the equations of 
elasticity, one numerical and the other semi-analytical. The latter method yields explicit 
formulas for stresses as expansions in Legendre polynomials about the solution for constant 
shell thickness. Though I only discuss the cases of a lithosphere thinner at the equator or at 
the poles, the method is applicable for any latitudinal variation of the lithospheric thickness. 
On lapetus, contraction or expansion on a lithosphere thinner at the equator explains the 
location and orientation of the equatorial ridge. On Mercury, the combination of contraction 
and despinning makes possible the existence of zonal provinces of thrust faults differing in 
orientation (north-south or east-west), which may be relevant to the orientation of lobate 
scarps. 
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1 Introduction 



A great variety of tectonic features is found on nearly all solid planets and large satellites of 
the Solar System: ridges and scarps, rifts and grabens, furrows and grooves etc. Their origin 
on Earth mainly lies in the movement of nearly rigid plates but other mechanisms must be 
found elsewhere since plate tectonics is unique to Earth. Though tectonic features can often 
be explained by a regional effect, such as an impact or the emplacement of a volcanic load, a 
subclass of them sometimes form a global pattern on the surface. In such a case, their cause must 
be the global deformation of the whole lithosphere which generates a global stress distribution 
resulting in a characteristic faulting pattern at the surface. 

Several processes affect the planetary figure. Though a change in the mean planetary radius 
(contraction or expansion) is the simplest conceivable deformation, its underlying causes can 
be complex since they d epend on the complica t ed ph ysics of th e interior of a self-gravit ating 
differentiated body (e.g. \Andrews-Hanna et al. 200^ for Mars, Sauvres and Croft 1986 | and 



Collins et al. 2009l | for icy satellites). On a planet with a lithosphere of constant thickness, con- 



traction and expansion bring about a ho mogeneous dist ribution of randomly o riented compres- 
siona l and extensional faults, respect ively Me/os/J . Il977l. The surfaces o f Mars [Anderson et al 



200l[ l2008t iKnavmever et all |2006| and Mercmv \Watters et all |2009| show widespread com- 
pressional tectonic features, termed wrinkle ridges and lobate scarps. Their distribution is 
however far from uniform and their orientation is not random. On Mars, regional effects 
such as the Tharsis rise have strongly influenced the global pattern of wrinkle ridges and 
the possible role of a glo bal contraction event remains under discussion Manaold et all 200d : 
Nahm and Schultz . 201Cl| . On Mercury, the lobate scarps that have been identified in Mariner 10 



images have a greater cumulative length in the southern latitudes, gener ally trend with i n 50 
of the north-south d i rectio n and preferably dip northward below 50° S 
Watters and Nimmo. l200{ 



Waiters et all 12004 : 



The anisotropy of lobate scarps is often attributed to the addi- 
tional contribution of des pinning (see below), but the regional effect of the Caloris basin has 



also been invoked [Thomas et all 1 19881] . A more complete and less illumination-biased fault 



catalog should result fro m the analysis o f Messenger images and alti metry data; results from the 



first flyby are promising Solomon et al . 2008t Watters et al . 2009 1. Among large icy sateUites, 
Ganymede is a showcase for extensional tectonics [Pavvalardo et al I [2004] . Global expansion 
is also though t to have played a r ole in the forination of the g lobal tectonic grids on Rhea 
IConsolm aand . 1985t Moore et a l.. 'l985; Thomas . 1988 1. Dione Moore . 1984 : Consolmaano . 
[l98 5] and Ariel Plescid . 19871 : \Nvf feneaaer et all 1997| . though the orientation of faults has 



been shown to be non-random in each case. Global contraction may a l so have occurr e d on R hea 



and Di one. New tectonic an alyses of Rhea \Moore and Schenk, '200?; ' Waaner et^alj, l2007l | and 
Dione Coff-Pochat and Collins. . ,20091: Waaner et al. . .2j30,9 : Stephan et al., 2010] are underway 
with Cassi ni data. While partial i magin g makes it difficult to identify a global tectonic grid 
on lapetus Singer and McKinnon . 20081. its huge equatorial ridge must be related to a global 



deformation. Contraction \Castillo-Roaez et al. . ^07] has been cited as a possible culprit, but 



the corresponding faulting pattern neither predicts the equatorial location nor the east-west 
orientation of the ridge. 

The next simplest deformation is a change in the planetary fiattening due to a decrease in 
the rotation rate, or despinning, that results from tidal effects leading to a resonant or syn- 



chronous rotation Murray and Dermott . 1999| . Despinning tectonics on a planet with a thin 



lithosphere of constant thickness consist of an equatorial zone of strike-slip faults (striking at 
abou t 60 degrees from the north) and of polar zones of extensional faults or joints striking east- 
west BurnA . Il976l : \Melosm . Il977f . D espun bodies are common in the Solar System: Mercury 
is in resonant rotation with the Sun Coldreich and Peale , Il968j l and nearly all large satellites 
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in the Solar System are in synchronous rotation with their parent body Peale . 1999| . Simul- 
taneous despinning and contraction h ave been used to justify the predominant north-sou t h ori- 
entation of lobate scarps on Mercury Melosh and Dzurisin , 1978t \Pechmann and Melosh , 1979t 



Dombard and Hauck\ . \200&^ . Beside the pattern of young lobate scar ps, M ercury exhibits a global 
grid o f more ancient lineaments which \ Melosh and DzurisiA [l978| and Pechmann and Melosh 



I979I interpreted as ev i dence of despinning (plus contraction), though the case is not closed 



Melosh and McKinnoA . Il988t \Thomas et all Il988| . While despinning is a generic phenomenon 



for large icy satellites, none exhibits unambiguous evidence of the corresponding global tectonic 
pattern, eithe r because of later resurfac ing or because despinning occurred before the surface was 
fully forme d ISauvres and Croftl. Il986i. The location of lapetus' ridge suggests a relation with 



despinning \Porco et all . 2005t Castillo-Roaez et al . 2007 1 but it was immediately noted that de- 
spinning tectonics cannot produce east-west features at the equator. Global tectonic patterns can 
also be produced by other types of deformations not examined i n this paper: reorientation rela - 
tive to th e spin axis \Melosh . 1980b : Leith and McKinnori. 1996 : Matsuvama and Nimmo . 200S 



Me/os/J. Il980at IffeZfews^em and Parmen^ieri.ll983l. diurna l tides and non-sync hronous 



Helfenstein and Parmentier . 19851 : Greenhera et al. . 1998t Wahr et al. . 2009| . These 



recession 
rotation 

deformations can be seen as a superposition of biaxial deformations and thus share the basics of 
the despinning model. 

Up to now the distribution and orientation of tectonic features have always been computed 
with the assumption of a lithosphere of constant thickness, one reason being the lack of data about 
the lithospheric thickness and its variations. The choice of a specific variation of lithospheric 
thickness may thus seem ad hoc if its cause is a variation in the internal heating flux. Two external 
phenomena provide a better motivation for a model of lithospheric thickness variation. On the 
one hand, the latitudinal variation in the solar insolation elevates the lithospheric isotherm at the 
equator of a planet without atmosphere. For Mercury, the lithosphere could be thinner by a factor 



of tw o at the equator in comparison with the poles [McKinnom . 119811 : \Melosh and McKinnon 



1988|; a longitude variation is also expected because of the difference in average temperature 



between the so-called equatorial hot and warm poles. On the other hand, tidal heating acting 
on the sat ellites with an eccen tric orbit creates additional lateral variations in the lithospheric 
thickness. 



Nimmo et al\ |2007j h owever found no evidence in the g lobal sha pe of Europa for th e 



thickness variations predicted by \Oiakanaas and Stevenson 1989l | (see also Tohie et al. 2003 1 ) 



Since large satellites are at present in a synchronous state of rotation, tidal heating is localized 
in longitude. Nevertheless tidal heating is evenly distributed in longitude before tidal locking 
is achieved, so that the lithospheric thickness only depends on latitude during the despinning 
phase. 

The goal of this paper is to predict global tectonic patterns due to contraction and despinning 
on a lithosphere with latitudinal thickness variation. Though I will assume in the numerical 
examples that the thickness variation is symmetrical about the equatorial plane, the method is 
also valid for non-symmetrical variations. The prediction of tectonic pattern s follows the fol l owing 



procedure. Fi rst the lithosphere is modeled as a thin elastic spherical shell [Turcotte et all 11981 



Beutha . l2008l | . Second the contraction or despinning is represented as a static deformation of 
the shell given by the difference between the initial and final figure of the planet. The stresses 
caused by the deformation of the thin elastic spherical shell are computed with the equations of 
equilibrium and elasticity. Third th e style and orientat ion of faults is predicted from the stresses 



with Anderson's theory of faulting Jaeaer et al. . 2007| 



The most significant result of this paper concerns the contraction of a planet with a lithosphere 
thinner at the equator (predictions for expansion follow by changing the sign of the stresses). 
If the contracting lithosphere is of constant thickness, stresses are homogeneous and isotropic, 
which means that their orientation and magnitude are independent of position and that the 
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principal tangential stresses are equal. The thinning of the lithosphere at the equator has two 
effects: (1) the stress becomes most compressive at the equator and (2) the meridional component 
of the stress becomes more compressive than the azimuthal component. This situation, as well 
as the corresponding expansion case, favors the development of faults striking east-west and 
starting preferably at the equator. The contraction of a shell thinner at the equator is thus 
complementary to despinning, for which the azimuthal stress is always more compressive than 
the meridional stress. If the inverse thickness is at most of harmonic degree two I show that this 
complementarity is described by an exact relation between a contracting planet and a planet that 
is spinning up. As for despinning, its tectonic pattern is not much affected by the variation of the 
lithospheric thickness: the main effect is the reduction in size of the polar provinces of tensional 
stress as the lithosphere thickens at the poles. The combination of contraction and despinning 
renders possible the existence of provinces of thrust faults differing in their orientation (north- 
south or east-west), but this only occurs for a specific ratio between contraction and despinning. 
If the lithosphere is thinner at the poles, global contraction (resp. expansion) leads to thrust 
(resp. normal) faults striking north-south and preferably formed in the polar regions, whereas 
the despinning pattern is weakly affected. 

As for applications, I show that contraction or expansion combined with lithospheric thinning 
at the equator provides an explanation for the location and orientation of lapetus' equatorial 
ridge. I also attempt to account for the orientation of lobate scarps on Mercury. Tectonic 
patterns resulting from simultaneous contraction and despinning on a lithosphere thinner at the 
equator are consistent with the latitudinal dependence of lobate scarps. The amount of required 
contraction is however much larger than what has been estimated from observations, making it 
necessary to resort to more contrived scenarios involving fault reactivation. 

The paper can be read at two levels: I present general results and geophysical interpretations 
in the main text, while the semi- analytical method of solution is described in detail in the 
appendices, which also include explicit solutions as well as a discussion of the relationship between 
thin shells, thick shells and Love numbers. Section[2]states the membrane equations that must be 
solved in order to compute the stresses generated by the deformation of a thin elastic shell with 
variable thickness. Remarkably, assumptions of axial and equatorial symmetry (presented in 
Section [3]) are sufficient for the prediction of the basic characteristics of stress and strain. These 
symmetries lead to dualities between the contraction and despinning solutions when the inverse 
thickness is at most of degree two, reducing by half the computational load and simplifying the 
discussion of combined contraction and despinning. I also compute a first order approximation of 
the contraction stresses from duality arguments. In Sectional I show how to solve the membrane 
equation. After discussing the parameterization of the thickness variation, I frame the problem 
for a numerical solution with Mathematica, before explaining a semi-analytical method in which 
truncated Legendre expansions serve to rewrite the membrane equation as a matrix equation. In 
Section [51 I use the contraction and despinning solutions of the membrane equation in order to 
predict the extent of tectonic provinces, i.e. areas with faults of given style and orientation. The 
cases of expansion and spinning-up are related to contraction and despinning by a sign change. 
In Section |6l I discuss the application of the predicted tectonic patterns to lapetus' ridge and 
Mercury's lobate scarps. 

2 Membrane stress in a deformed spherical shell 

The style and orientation of faults due to contraction and despinning can be predicted from 
knowledge of the stresses in the lithosphere. The lithosphere is modeled as a thin elastic spherical 
shell (in this paper, lithospheric thickness and elastic thickness are synonymous). The important 
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assumptions that the sheh is thin and that it is elastic are analyzed h\ \Melosh 1977 1 for a shell 
of constant thickness. 

First, Melosh shows that models with a thin or a thick shell predict very similar despinning 
faulting patterns, except the possible presence of an equatorial thrust fault province when the 
shell is thick (the contraction tectonic pattern is not affected by the thickness of the shell). 
Despinning stresses at the surface of a thick shell can easily be computed if the secular tidal Love 
numbers ft.^ and I2 are known (see Appendix 18. 4p . At a given latitude, the stress magnitude 
mainly depends on /i^ whereas the relative size of the principal stresses depends on the ratio 
^2/^2- When the shell thickness increases, /i^" decreases so that despinning stresses decrease 
in magnit ude. While Melosh computes t hin shell stresses under the assumption of hydrostatic 
flattening, Matsuvama and NimmJl 2008j extend the domain of application of thin shell formulas 
by using a non- hydrostatic flattening parameterized by /i^". This latter procedure is equivalent 
to keeping the overall /i^-dependence in thick shell formulas while setting the ratio I2 /h'^ to its 
membrane value (1 + v)/{J> + v), where v is Poisson's ratio (see Appendix 18. 4p . 

Second, Melosh argues that the faulting pattern predicted by the elastic model is not sub- 
stantially modified when the lithosphere becomes plastic, though the numerical values of the 
stresses are bounded by the yield stress. The elastic model is thus limited to the prediction of 
the faulting style (including its orientation) but does not say whether faulting occurs. 

Stresses in a shell satisfy equilibrium equations which relate them to the load deforming the 
shell. However the load causing the contraction or the flattening change during despinning is a 
priori unknown. Thus the full formulation of elastic equations involving stress-strain and strain- 
displacement relations becomes nece ssary. For a n elastic shell of constant thickness, contraction 
stresses were obtained by Lame [see Love . 1944 p. 142] whereas despinning stresses were com- 
puted by I Venina-Meinesz 1947 1 . There are no analytic fo rmulas giving equivalent results for a 
thin elastic shell with variable thickness but BeuthA |2008l | recently derived the scalar equations 
governing the deformation of such a shell under an arbitrary load. 

Since displacements of the shell surface due to contraction and despinning have a large wave- 
length, bending moments in the shell are negligible: the shell mainly re sponds by stretching and 



198lj . In the membrane limit. 



is said to be in the membrane regime \Kraus , 1967 ; Turcotte et al. 
bending moments are set to zero: deformation equations take a simplified form which can be 
obtained from the full deformation equations by setting to zero the bending rigidity D. In that 
case, the membrane equations for a spherical shell of radius R and variable thickness h are 



A'F 

C{a;F)-{l + :y)A{a;F) 



-Rq, 



(1) 
(2) 



where q is the transverse load (positive inward), w is the transverse displacement (positive 
outward) and F is an auxiliary function called the stress function. Tensile stress is positive. 
The tangential load has been set to zero. The differential operators A', C and A are defined in 
Appendix 18.1 1 bv Eqs. (|63I) - ([S7)) . The thickness h is included in the elastic parameter a: 



1 

Eh 



(3) 



E is Young's modulus (which is allowed to be variable), v is Poisson's ratio (which must be 
constant) . 

In the membrane limit, the stress integrated over the thickness of the shell, or stress resultant, 
is given by 

iNgg,N^^,Ng^)^i02,Oi,-03)F. (4) 
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with the operators Oi being defined by Eqs. The stress as expressed in Eq. 

automaticahy satisfies the two tangential equihbrium equations while the radial equilibrium 
equation is encoded in Eq. ([T]). In the membrane limit, there is thus only one equation where 
elastic parameters appear, Eq. ([2]), which results from an equation of compatibility between the 
strains. 

In the membrane limit, the stress is constant through the shell thickness. It is thus related 
to the stress resultant by 

c^-y = j^N.j . (5) 



The strain is given by 



1 

E 
1 



= [cree - , (6) 



^ifip = (o-ifip - i^cree) , (7) 

If the load is known, the stress function F can be computed from the first membrane equation 
([T]) and the stress from Eqs. (l4])-([5]). In that case, the elastic thickness only appears in Eq. ([5]), 
which means that the solution for variable thickness is directly obtained from the solution of 
constant thickness by the scaling ho/h (ho is the mean elastic thickness). This procedure implies 
that the load is the same whether the elastic thickness is constant or not. However the load 
is a priori unknown when modeling contracting and despinning events; instead, the transverse 
displacement w is the known input. The stress function F must thus be computed from the second 
membrane equation, Eq. ([2]), which is linear in F and w with variable coefficients depending on 
1/ Eh. As a corollary, the load causing contraction or despinning is not the same in the cases of 
constant and variable elastic thickness. 

It should be noted that the variation of the lithospheric thickness has an effect on the de- 
formation. Contraction is not a problem if it is defined as a uniform change in radius, but it is 
another matter if the definition involves physical processes such as an internal thermal contrac- 
tion or an expansion due to water-ice transition. Despinning causes deformations of harmonic 
degrees zero and two (the former is very small and usually neglected) if the lithospheric thickness 
is constant but other harmonic degrees appear if it is variable. In principle the effect of internal 
contraction and despinning on the shape of a planet with a variable lithospheric thickness should 
be computed with a three-dimensional model of the interior. Love numbers quantify the response 
of a planet with a spherically symmetric internal structure (see Appendix 18. 4p , but there is no 
standard technique dealing with the non-spherical case. Extreme cases of a very thin or a very 
thick lithosphere are not a problem. On the one hand, the lithosphere does not affect the shape 
of the body in the limit of vanishing thickness. On the other hand, thickness variations of a thick 
lithosphere are expected to be relatively small with respect to the average thickness: latitudinal 
variations in solar insolation, for example, induce smaller thickness variations (in percentage) in 
a thick lithosphere than in a thin one. In both cases, the response of the planet is thus adequately 
described by using a model with a spherically symmetric internal structure. Between these ex- 
tremes, the size of the deformations with harmonic degrees not equal to zero (for contraction) or 
two (for despinning) should be investigated, but this question is beyond the scope of this paper. 
I thus choose to define contraction and despinning as phenomena associated with deformations 
of harmonic degrees zero and two, respectively. 

How good an approximation is the membrane limit? If bending terms are not neglected, 
Eq. ^ remains the same but new terms appear in Eq. Q for the stress resultants. These new 
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terms have the form {D/B?) OiW (with D = Eh^ /12{1 — v"^)) so that they are smaher than the 
dominant term by a factor ^ = 12B?'/h?. Neglected bending terms are thus of the same 

order of magnitud e as oth e r sma U terms present in the equations before the 'thin shell' limit 
of large ^ is taken [Beuth^ . |2008 |. It is possible to solve the complete equations (with bending 



terms and no large ^ limit) with the same method as proposed in this paper, but the gain in 
precision would only be of order 1/^. Even so the theory remains a thin shell theory in the sense 
that the radial stress is neglected. 

The membrane and stress equations ©-© become nondimensional with the following defi- 
nitions: 



a 


a 

ao ' 


(9) 


F 


= o^oF , 


(10) 


w 


w 
R ' 


(11) 






(12) 




1 


(13) 



where ao is the average of a on the sphere, i.e. the coefficient of degree zero in the expansion of 
a in Legendre polynomials (ao 7^ "571^ )■ ^"^^ membrane equation ([2]) then becomes 

C{a;F)-{l + iy)A{a;F) ^A'w. (14) 

The nondimensional stress is computed with the nondimensional version of Eqs. Q-®: 

{Ng0,N^^,Ne^) = (02,Oi,-03)i^, (15) 
CTy = aNij. (16) 

The strain is already nondimensional: 

eee = a{02-iyOi)F, (17) 
e^^ = a{0i-v02)F, (18) 
eg^ = -a{l + i^)03F. (19) 

The solution of the nondimensional membrane equation (jl4p when w is only of degree zero 
{w = wqPq) will be called the contraction solution F^ (whatever the sign of wq), whereas the 
solution when H) is only of degree two {w = 'W2P2) will be called the despinning solution F^ 
(whatever the sign oi W2)- Negative (resp. positive) values of wq correspond to contracting (resp. 
expanding) planets, whereas positive (resp. negative) values oiw2 correspond to despinning (resp. 
spinning-up) planets. Solutions for a constant elastic thickness are computed in Appendix [ 



3 Symmetry and duality 

As discussed in the introduction, I make the simplifying assumption that the value of the elastic 
thickness only depends on the latitude. The contraction and despinning problems are then 
axially symmetric. Though the numerical method developed in Section S] does not require it, I 
also assume equatorial symmetry, i.e. mirror symmetry about the equatorial plane. These two 
assumptions approximately hold for the physical causes of the thickness variation considered here: 
variation in solar insolation and tidal heating during despinning, provided that the obliquity of 
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(a) EXPANSION (b) DESPINNING 




-1 I I I I I i-J — u.uuc-i I 1 1 1 1 lj 

15 30 45 60 75 90 15 30 45 60 75 90 



Colatitude IdegreesJ Colatitude [degrees] 

Fiff Ur6 li Strain as predicted from symmetry arguments: (a) expansion, (b) despinning. Continuous curves refer 
to egg while dashed curves refer to £ipip. Three cases are shown: constant thickness, equatorial thinning, polar 
thinning. The unlabeled middle curves in the despinning case are for constant elastic thickness. The numerical 
values used to draw the curves are: tDo = 0.001 for expansion, tD2 = 0.001 for despinning, r = 1/2 for equatorial 
thinning, r = 2 for polar thinning (r is defined by Eq. II34II ). Poisson's ratio equal to 0.25. Extensional strain is 
positive. 



the body with respect to the Sun (in the former case) or to the tidal source (in the latter case) 
is small. 

Using symmetry and duality arguments, I describe in this section the characteristics of stress 
and strain on contracting or despinning shells thinner (or thicker) at the equator. The duality 
relation will also serve to obtain a numerical approximation of the contraction solution. 

3.1 Axial and equatorial symmetry 

The main characteristics of the contraction solution can be determined from symmetry argu- 
ments, without solving the membrane equation. Axial symmetry gives the following constraints 
on stress and strain (see also Fig.[T|): 

1. the meridional and azimuthal stresses (resp. strains) are principal stresses (resp. strains), 
i.e. asip = (resp. eg^, = 0). 

2. the meridional and azimuthal stresses are equal at the poles; this is also true of the strains. 

3. the slopes of the stress and strain components vanish at the poles (it is a consequence of 
the previous point combined with the equations of equilibrium); they also vanish at the 
equator if there is equatorial symmetry. 

4. the average of the meridional strain is independent of the elastic thickness (see Eq. (IllBp ): 
< egg >= wq for contraction, < egg >— ■iZ'2/4 for despinning. 

5. the azimuthal strain at the equator is independent of the elastic thickness (see Eq. (jllSp ): 
^ipip\e=f = Wo for contraction, e^p^plg^:^ — —W2/2 for despinning. 

These constraints allow me to predict the response of the shell to a radius variation. The 
stresses for a constant elastic thickness are given in Appendix 18.31 Fig. [T] illustrates the following 
discussion. I suppose that the shell is expanding (wq > 0) so that stresses and strains are positive 
(the solution for a contracting shell is obtained by changing the sign of stresses and strains). I 
also assume equatorial symmetry with a monotonous variation of the thickness from the pole to 
the equator, so that stress and strain components can be expected to vary monotonously. Stress 
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and strain are concentrated where the shell is weaker. In other words, the membrane stretching 
is maximum where the shell is thinnest. If the shell is thinner at the equator, the azimuthal 
strain is thus maximum at the equator, where it takes the value wq imposed by the equatorial 
length variation, and minimum at the poles, where it is equal to the meridional strain. The 
meridional strain is also maximum at the equator and minimum at the poles, taking somewhere 
in between its average value wq. Therefore, the meridional strain always exceeds the azimuthal 
strain when the shell is thinner at the equator. This ordering is conserved for the stresses: the 
meridional stress always exceeds the azimuthal stress when the shell is thinner at the equator, as 
shown by the inversion of Eqs. ([6|)-([8]) with v < I. If the shell is thinner at the poles, a similar 
analysis leads to the conclusion that the strain (resp. stress) is maximum at the poles and that 
the azimuthal strain (resp. stress) always exceeds the meridional strain (resp. stress). 

The despinning shell case {iju2 > 0) is not as simple since the combined effects of the non- 
spherical deformation and thickness variation can lead to local extrema of the stress and strain 
components between the poles and the equator. However the ordering of the strains (e,^,^ more 
compressive than egg) is expected to be the same as in case of constant elastic thickness, since 
the meridional strain has an average positive value (w;2/4) while the azimuthal strain varies 
between the value of the meridional strain at the poles and the negative value — u'2/2 at the 
equator. The location of the shell thinning has thus less influence on the stress/strain pattern 
than in the contraction case. Because of stress/strain concentration in weaker areas, the stress 
and strain components are expected to come closer to (resp. farther from) zero at the poles 
if the shell becomes thinner at the equator (resp. at the poles), as shown in Fig. [TJd. The 
meridional strain curve will accordingly move so as to keep its average constant. The principle of 
stress/strain concentration does not apply to the azimuthal strain at the equator since its value 
is fixed by the geometry. As above, the ordering is conserved when going from strains to stresses: 
atptp is more compressive than agg. Another argument for ordering will be obtained from the 
contraction-despinning duality in Section \3l2\ 



3.2 Contraction-despinning duality 

Still assuming axial and equatorial symmetry, I now show that the contraction and despinning 
solutions are not independent if the harmonic expansion of the inverse thickness is limited to 
degree two. Even if this assumption is not strictly valid, the harmonic of degree two will dom- 
inate if the thickness variation is monotonous between the poles and the equator. Under these 
assumptions, a can be written in terms of Legendre polynomials of degrees zero and two: 

a = l + a2P2. (20) 

where P2 is the Legendre polynomial of degree two defined by Eq. (j7ip . The parameter 012 is in 
the range [—1,2] and is negative (resp. positive) for shells thinner at the equator (resp. poles). 
The parameterization of the thickness variation is discussed in more detail in Section 14.11 
The stress functions for contraction and despinning then satisfy the following duality: 

— = "2 — ■ (21) 

Wq \ — V W2 

The first term in the right-hand side is the contraction solution when the elastic thickness is 
constant. The duality (PT|) can be checked by substitution in Eq. ((TH) and using the identity 
Since spherical harmonics of degree one belong to the kernel of A' (see Appendix 18. II) . the 
duality remains valid if a includes such harmonics. If a contains harmonics up to degree fc > 2, a 
similar relation exists (see Eq. (1125^ ) between the solutions for right-hand sides of higher degree 
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(w = W3 ■■■ , w = Wk Pk), ^ use it in Appendix 18.61 for the computation of the contraction 
solution at first order. I also give below a relation of fourth degree between the stresses in order 
to discuss the combination of contraction and despinning. 

Eq. (l2lT) leads to a duality between the stress resultants for the contraction and despinning 
solutions: 

Wq 1 — 1/ W2 

The stress resultants for a contracting shell are thus the same as those for a despinning shell 
with W2 = —a2WQ, plus the stress resultants for a contracting shell of constant thickness. If the 
contracting shell (wq < 0) is thinner at the equator (a2 < 0), the equivalent 'despinning' shell 
is actually spinning faster since W2 < 0. This duality is clearly visible by comparing Figs. 
and |3)d. a similar duality holds for the stresses, but the first term in the right-hand side now 
depends on latitude: 

^ - a2 ^ . (23) 
Wo 1 — ly W2 

Finally the duality relation for the strains reads 

-^ = a6,,-a2^. (24) 

Wo W2 

These dualities are extremely useful. Their most obvious use is to generate both contraction 
and despinning solutions by computing only one of them, but they also have other advantages. 
First, Eq. shows that the ordering in size of the components of the stress is either the same 
(if a2 > 0) or reversed (if ck2 < 0) when going from the contraction to the despinning solution. 
Thus the azimuthal stress is more compressive than the meridional stress for a despinning shell 
whether the shell is thinner at the equator or at the poles. While the isotropic term Q.Sij/{l — v) 
in Eq. ([23| does not affect the ordering of the stress components, its spatial dependence influences 
the positions of the maxima which must thus be computed from the numerical solution. 

Second, the dualities serve to compute a first approximation of the contraction solution on 
a shell of variable thickness, knowing the despinning solution on a shell of constant thickness. 
The trick comes from the fact that the term linear in 612 in the contraction solution is related 
by duality to the term independent of 0,2 in the despinning solution. The contraction solution 
is computed in such a way in Appendix 18.61 (see Eqs. (|121l) - (|124p ). In dimensional notation, the 
contraction stresses at first order in 0,2 read 

c ^ _ /^T , o ^ + (2 + J^) cos 26* 



^ 2(5 + .) j' 

c ^ A^o- l + (l + 2.) cos 20 ^ 

where (Tq is the stress for a radial contraction of {SR > 0) when the thickness is constant: 

E SR 

Fig. 113b shows that the first order approximation is good if 0.2 — —2/5 and rather bad if 02 
is close to its lower bound —1. In Appendix 18.61 the method is generalized to an a having an 
arbitrary Legendre expansion (see Eqs. (|128l) - (|129l) '). Though the basic characteristics of stress 
and strain for a despinning shell could be predicted by the symmetry arguments of Section 13.11 
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duality cannot serve for the computation of their numerical values. Their formulas at first order 
in Oil are given in Appendix 18.101 

A third interesting consequence of the dualities concerns the analysis of combined contraction 
and despinning. If the amounts of contraction and despinning are related by 

W2 = a2Wo , (28) 

the total stress computed with the duality (|23l) is isotropic (though not homogeneous): 

a%^af^^aw^j^. (29) 

Consider a contracting (wq < 0) and despinning (w2 > 0) body with a lithosphere thinner at 
the equator (a2 < 0). If there is less (resp. more) despinning than the threshold the 
meridional stress is more (resp. less) compressive than the azimuthal stress. This fact is useful 
when determining the orientation of faults for contraction plus despinning. The threshold is 
reached for a contraction 6R {6R > 0) of 

SR=^ -J 1^ — , (30) 

Q!2 3g 

where Hi (resp. f2^) is the initial (resp. final) angular rate, g is the surface gravity and ft.^ is 
the degree-two displacement Love number. The above equation results from the substitution of 
Eqs. dHll) and ^ into Eq. 

The threshold (j28p becomes dependent on latitude when the inverse thickness expansion 
includes higher harmonics. If a = 1 + 012 P2 + 0^4 -P47 the generalization of the stress duality reads 

^+a2^+a4^=a-^, (31) 

Wo W2 W4 1 ~ V 



where the index E denotes the solution of the membrane equation (1141 for w = WiPi, which 
is given at zeroth order in (012 ,014) by Eqs. ([M)) - (P5|) . Since there is actually no deformation 
of degree four, the third term in the left-hand side of Eq. (I5T]) should not be interpreted as a 
physical stress but rather as a deviation from the duality At the threshold (PSI) . the total 

stress afj = af^ -f (j^^ is thus not isotropic anymore: 

i^^P^ 10^2 - 2) . (32) 

If d4^wo > 0, the meridional stress is more (resp. less) compressive than the azimuthal stress for 
high (resp. low) latitudes, with the turning point being given by a latitude of arcsin(l/\/7) = 
22.2° (this value is only valid if Q;2,a4 <C 1). If a^wo < 0, the high and low latitudes zones 
are exchanged. In Section I propose a parameterization of the inverse thickness which has an 
unlimited expansion in Legendre polynomials (see Eq. p6p ). Fig. IH] then shows the latitudinal 
dependence of the threshold between north-south and east-west thrust faults for different types 
of thickness variations. The case described by Eq. ([5^ corresponds to a value of the parameter 
k close to zero. 



4 Solving the membrane equation 

Symmetry and duality arguments have served us well for the determination of the main char- 
acteristics of the stress and strain curves, as well as for a first numerical approximation of the 
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contraction solution. It remains however necessary to solve the full membrane equation in order 
to compute both the contraction and despinning solutions at an arbitrary degree of precision. I 
will present two methods which are valid for an arbitrary deformation of the surface. The first 
one is numerical and has the advantage of minimizing the amount of programming. The second 
one is semi-analytical, in the sense that the solution can be expressed as a perturbation expansion 
in Legendre polynomials about the solution for constant shell thickness. It has the advantage 
of producing explicit formulas for the stresses in which the influence of parameters describing 
thickness variations clearly appears. Before tackling the methods of resolution, I briefly discuss 
the parameterization of the thickness variation. 



4.1 Thickness variation 

Since my aim is to predict tectonic patterns without assuming a specific lithospheric structure, 
I will work with a simple parameterization of the variation of the shell thickness. I choose to 
parameterize a instead of the thickness h because it is more convenient for the semi-analytical 
method. A few assumptions constrain the possible form of the thickness variation. Variations 
caused by solar insolation and tidal heating are of very large wavelength, so that only slow- 
varying functions of the colatitude are permissible. Axial symmetry means that the thickness 
only depends on the colatitude 9, but it also imposes that the derivative of the thickness vanishes 
at the poles. If there is equatorial symmetry, the derivative of the thickness vanishes at the 
equator (this restriction can be lifted without changing the methods of resolution). A simple 
parameterization of the thickness compatible with these constraints is given by Eq. (j20p : 

a(6l) = 1 + 02^2 (cos , (33) 

where the dependence on the colatitude 9 has been made explicit. I define the equator-to-pole 
thickness ratio r by 

r-^, (34) 
lip 

where He is the equatorial thickness and hp the polar thickness. The coefficient a2 is then given 
by 

a2=2^. (35) 

The thickness is positive everywhere if —1 < a2 < 2. Negative and positive values of a2 
describe elastic shells that are thinner at the equator and at the poles, respectively. Fig. [2] 
shows various profiles of a and h/hp when the lithosphere is thinner at the equator: the values 
r = (1, 1/2, 1/4, 1/10) correspond to a2 = (0, -2/5, -2/3, -6/7). 

Although the parameterization is sufficient for the determination of the most important 
features of tectonic patterns, it suffers from a weakness revealed by Fig. ^jp. When the equator- 
to-pole thickness ratio decreases, the zone of equatorial thinning expands toward the pole. In 
other words, the parameter a2 not only affects the equator-to-pole thickness ratio, but it also 
modifies the relative extension of the thin and thick zones. Therefore, the parameterization p3p 
is not appropriate for the analysis of the influence of the thin zone size. I thus deflne a new 
parameterization in which the interval [0, 7r/2] is non-linearly stretched by a function ipk- As a 
result, the extension of the thin zone can be modulated without affecting the equator-to-pole 
thickness ratio: 

a(0) = 1 + ^2^2 (cosVfc(^)), (36) 

where —1 < 52 < 2. The 'tilde' notation indicates that the normalization differs from the one 
used for a because the Legendre coefficient of degree zero differs from one. The term P2(cos ipk{S)) 
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(a) INVERSE THICKNESS (b) THICKNESS 




15 30 45 60 75 90 15 30 45 60 75 90 



Golatitude [degrees] Colatitude [degrees] 

Figure 2; Equatorial thinning of the Hthosphere as parameterized by Eqs. II33I I-I I35I I: (a) profiles of a (pro- 
portional to the inverse thickness), (b) thickness profiles normalized by the polar thickness. The four profiles 
correspond to equator-to-pole thickness ratios r equal to (1, 1/2, 1/4, 1/10). 



has generally an unlimited expansion in Legendre polynomials, including a term of degree zero. 
Nonetheless the normalization of a does not matter for the computation of stress and strain; it 
only affects the stress resultants. 

The function ipk{0) is defined on the interval [0, 7r/2] by 



7r sinti(2fc9/7r) 



if fc > 



M0) = \ I tanL-^i^^/.) .f^^o ' (37) 
2 ta,nh(k) 

It can extended to the interval [0,7r] by mirror symmetry about the equatorial plane. The limit 
fc "> yields the parameterization ([55]) : 

limiPkiO) =0- (38) 

fe->0 

Fig. [3] shows that positive values of k correspond to a reduction in the size of the thin zone. The 
equator-to-pole thickness ratio is left unchanged since it is controlled by a2- Conversely, negative 
values of fc correspond to an extension of the thin zone. 

4.2 Numerical method 



The software Mathematica [Wolfram Researcm . l2008l | has a powerful command called 'NDSolve' 



for the numerical solution of differential equations. When there is axial symmetry, the second 
membrane equation (1141) is an ordinary differential equation of order four. I found it convenient 
to solve it as a system of two differential equations of order two: 

A'F = -q, (39) 
-A' {aq) - {l + v)A{a]F) = A'w . (40) 

where the nondimensional load q is defined by 

q = RaQq. (41) 

The various terms in the above equations can be expanded with the formulas of Appendix 18.11 

A'F = AF + 2F, (42) 
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(a) INVERSE THICKNESS 



(b) THICKNESS 




15 30 45 60 75 90 15 30 45 60 75 90 

Golatitude [degrees] Colatitude [degrees] 



Figure 3: Equatorial thinning of the hthosphere as parameterized by Eqs. Il36l l-l l551l : (a) profiles of a (propor- 
tional to the inverse thickness), (b) thickness profiles normalized by the polar thickness. The equator-to-pole 
thickness ratio r is equal to 1/2: the curves fc = are the same as the curves r = 1/2 in Fig. [2] Curves A; = 2,4 
correspond to a reduction of the thin zone whereas curves k = —2, —4 correspond to an extension of the thin 
zone. 



A'(ag) = (Aa)g + 2— ^+a(Ag) + 2a9, (43) 

with 

In principle the system ([Mjl - lHUl) should be solved on the segment [0, tt] but the points and 
TT must be excluded from the interval because of the apparent polar singularities in spherical 
coordinates. I thus restrict the interval to [a, b] with a = e and 6 = tt — e (e is a small positive 
number, for example 0.001). 

Four boundary conditions are required. Smooth solutions with axial symmetry satisfy 

-QQ\e=a = ^|e=a=0, (46) 

Before applying these equations as boundary conditions, I should make a caveat. The stress 
function has the particularity that it can be redefined at will because of the degree-one freedom 
mentioned in Appendix 18. II 

F F + PcosB , (48) 

where /3 is an arbitrary real number. Since cos 6* satisfies Eqs. (|15)) - (H7)) in the limit e — > 0, 
another condition must be specified in order to fully determine F, for example: 

F|e=a = 0. (49) 

This last condition being arbitrary, one should keep in mind that the same physical problem can 
be solved with various boundary conditions. The resulting stress functions only differ by their 
degree-one content, i.e. by cos 6* multiplied by some constant. The corresponding stresses are of 
course identical. 
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Suitable boundary conditions for the interval [a,b] consist of Eq. (HH]), plus three conditions 
chosen among Eqs. P5)) - (|T7)) . I found it best to specify three conditions on F and one condition 
on q. In most cases NDSolve also works well with only Eqs. (P5| - (j47p as boundary conditions, 
meaning that its algorithm chooses one possible solution for F among an infinity (note that the 
freedom in the definition of F only appears in the limit e — > 0). Though it is not necessary, one 
can get rid of the cos 9 component of F by orthogonalization: 

F^F-\[ Fcosesinede. (50) 

If there is equatorial symmetry, it is sufficient to consider the interval [a, h] with a = e and 
6 = 7r/2. In that case, Eqs. (j46p - (|47)) form a suitable set of boundary conditions since cos 61 
violates Eq. (|T7)l when b = 7r/2. The cosO component of F then vanishes. 



4.3 Semi-analytical method 

If the elastic thickness is constant, the membrane equation ([H]) is diagonal in the basis of spherical 
harmonics. The contraction solution is then of harmonic degree zero whereas the despinning 
solution is of degree two (see Appendix 18. 3|) . This straightforward method cannot be used when 
the elastic thickness is spatially variable because of the coupling of the different harmonic degrees 
and orders. Here I show that the system of coupled differential equations remains manageable 
under the assumption of axial symmetry. Since the imposed deformation and the elastic thickness 
do not depend on the longitude, only zonal spherical harmonics contribute. 

The problem thus consists in solving the membrane equation (jl4p with expansions in zonal 
spherical harmonics, i.e. Legendre polynomials: a = ^ap P^, F = '^FgPg and w = '^wgPi. 
The action of the operators C and A on Legendre polynomials produces a finite sum of Leg- 
endre polynomials which is computed in Appendices 18.71 and 18.81 (final results are embodied in 
Eqs. (|146|) - (|147p ). If the expansions are truncated at degree n, the membrane equation takes a 
matrix form with each row corresponding to a harmonic degree on which the membrane equation 
is projected, except for degree one: 

]V[(")f = A'w, (51) 
with the membrane matrix M^"-' being defined by 

M(") = C(") - (1 + I/) A(") . (52) 

C*^") and A*^"' are matrices approximating at degree n the operators C and A, respectively. F 
and w are vectors containing the coefficients Fi and wi, respectively. A' is the diagonal matrix 
with elements = —£{£ + 1) + 2. The absence of a row for degree one is due to the fact that 
spherical surface harmonics of degree one do not belong to the images of the operators C and A 
(see Appendix 18. 8p . Correspondingly, there is no coefficient of degree one in the vectors F and w 
since the membrane equation does not constrain them (see Appendix 18. ip . Given a deformation 
w, the membrane equation in its matrix form can for example be solved by matrix inversion, 
yielding the Legendre coefficients of the stress function: 

F = (m("))"^ A'w, (53) 

An example of the membrane matrix and of its solution is given in Appendix 18 . 91 for an expansion 
of a limited to degree two and a truncation degree equal to 6. If the thickness variation is 
symmetric about the equatorial plane, only Legendre polynomials of even degree will contribute 
to the solution because the deformation is also symmetric (degree zero or two). 
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Since a = 1 + 02^2 is finite even if hs or hp vanishes (though not both), the nondimensional 
membrane equation can in principle be solved with a vanishing thickness at the equator or at 
the pole. Problems of divergence however appear in these extreme cases for the nondimensional 
stress function and the stress resultants, though the stresses and the strains are well behaved. 
The membrane matrix M*^"^ is invertible in the physical range of v € [0, 1/2] and for 0:2 & [—1, 2], 
as can be seen from the examination of its eigenvalues. 

If a is limited to degree two, the approximate membrane equation (|5ip can also be solved by 
expanding F in a2 and solving order by order: 

00 

F = ^(a2)PF(P). (54) 

p=0 

It is convenient to split the membrane matrix so that the dependence in 02 becomes explicit: 

M(") ^ M*,") + 02 M^"^ . (55) 

(n) 

The matrix Mq is diagonal so that its inversion is straightforward. The solution is then gener- 
ated order by order in 0:2 by a recurrence relation: 

F(p+i) = (m^"^) M^"^ F(p) , (56) 

which is initiated with F*^°\ the solution of Mq"''f = A'w. This method can be generalized to 
an a of degree higher than two. Since the largest eigenvalues of the iteration matrix appearing 
in Eq. (155)) tend to 1 (from below) as the truncation degree increases, the series ([5^ may diverge 
for near extremal values of ck2. If 1 0:2! < 1, the series converges but the convergence can be very 
slow if 0:2 is close to -1. Even if it is not the best numerical method, the cii-expansion provides an 
explicit solution which is helpful to understand properties such as the rule governing the decrease 
of the Legendre coefficients of the solution (see Appendix 18. 9|) or the pseudo- nodes in the stress 
and strain components (see Appendices 18.61 and I8.10p . 



5 Tectonic patterns 



5.1 Stress and faulting 

Tectonic patterns can be predicted from the analysis of stresses. According to Anderson's theory 
of faulting, the faulting style dep ends on how the vertical (or radial) stress conipares with the 
horizontal (or tangential) stresses [ Melosh . 1977t [jaeqer et al. . 2007 : Schultz et ai . 2009l |. Thrust 
faults, strike-slip faults and normal faults occur if the radial stress is respectively least compres- 
sive, intermediate or most compressive among the principal stresses. Thrust faults and normal 
faults strike in the direction of the intermediate principal stress, whereas strike-slip faults strike 
in a direction at about 30° from the direction of the most compressive s tress. The radial stress 

An- 



vanishes since faults occur at the surface (this assumption is criticized hv \GolombeM 1985 
derson's theory presupposes that all principal stresses are compressive, at least when including 
the lithostatic pressure that must be added to the stresses computed from thin shell theory. I 
will briefly mention below the possible occurrence of near-surface tensile fai lure when one (or 
more) principal stress is extensional. This point is discussed in more depth hY lMelosh 1977 [ and 



Schultz and Zuhen |1994 | 



Once the contraction solution of the membrane equation ^ has been found with the methods 
of Section it is possible to compute stress resultants, stresses and strains with the formulas 
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Figure 4: Stress resultant as a function of colatitude when the lithosphere is thinner at the equator: (a) 
contraction of 0.1% (iDo = —0.001), (b) despinning with flattening change 5f = —0.0015 (■102 = 0.001). Continuous 
curves refer to Ngg while dashed curves refer to Ntp^. Elastic parameters are E = 100 GPa and v = 0.25. The 
thickness variation is parameterized by the equator-to-pole thickness ratio r as in Fig. [21 the average inverse 
thickness is equal to 1/(100 km). Truncation degree n is equal to 20. Tensile stress is positive. The duality II22II 
is apparent when comparing contraction and despinning curves. 



of Section [2l The nondimensional stress is sufficient for the determination of tectonic patterns 
since the physical magnitude of the stress is only useful when comparing the predicted stresses 
with a failure criterion. Nonetheless figures represent dimensional quantities for values of param- 
eters typical of terrestrial planets: E = 100 GPa, v = 0.25, average inverse thickness equal to 
1/(100 km), contraction of 0.1% and flattening reduction of 0.15%. I initially give examples with 
a limited to harmonic degree two, i.e. a is parameterized by Eq. p3l) . Most basic characteristics 
of the tectonic patterns are indeed independent of the precise latitudinal variation of the elastic 
thickness as long as it smoothly varies from the pole to the equator. Figs. 2] and O show the 
stress resultants and the stresses for the contraction and despinning of a lithosphere thinner at 
the equator. The duality between contraction and despinning solutions appears clearly on Fig. H) 
The stress curves for varying a2 seem to cross at a fixed point on Fig. [5] but this is only true at 
first order in a2 (see Appendices 18 . 61 and 18 . lOp . The case of polar thinning is illustrated by Fig.[Tl 
which shows the strains for expansion and despinning (stress curves are very similar). Besides I 
will resort to the parameterization p6p in order to demonstrate two effects due to the variation 
in size of the thin zone: it modifies the position of the extremum of the meridional stress and 
it affects the orientation of thrust faults when contraction is combined with despinning. Fig. [6] 
shows the stresses for contraction and despinning when the thin zone size is extended or reduced. 

5.2 Contraction or expansion 

For the contraction of a shell thinner at the equator, the stresses have the following characteristics: 

• Stresses are everywhere compressive. 

• The meridional stress is more compressive than the azimuthal stress (Fig. [5^); This situa- 
tion favors the development of thrust faults striking east-west. 

• The stress is most compressive at the equator (Fig. [5^). This situation favors the develop- 
ment of thrust faults at the equator. Note that the stress resultant is most compressive at 
the poles (Fig. Hi). 

The contraction fault pattern predicted with Anderson's theory is thus significantly modified by 
a variable elastic thickness: instead of thrust faults with random strikes, the predicted pattern 
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Figure 5: Stress as a function of colatitude when the Uthosphere is thinner at the equator: (a) contraction, (b) 
despinning. Continuous curves refer to agg while dashed curves refer to a^^. Other details as in Fig.U] 



(a) CONTRACTION (b) DESPINNING 




15 30 45 60 75 90 15 30 45 60 75 90 

Colatitude IdegreesI Colatitude Idegrees] 



Figure 6: Influence of thin zone size on stress: (a) contraction, (b) despinning. Continuous curves refer to agg 
while dashed curves refer to cr^^. The thin zone size varies as in Eq. II36I I (see also Fig. [3] but r is different). The 
equator-to-pole thickness ratio r is equal to 1/10 {a2 = ~6/7). The curve k = 2 (resp. k = —2) corresponds to 
a thin zone reduction (resp. enlargement). The curve fe = is the same as the curve r = 1/10 in Fig. [5] Other 
details as in Fig. |4] 
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consists of thrust faults striking east-west and preferably formed at the equator. The third 
characteristic is not strictly true for all thickness profiles. The curve k = —2 in Fig. [S^ shows 
that two things happen if the thin zone has a large extension. First stress curves become nearly 
flat in that zone, where the stresses approach the value predicted for a lithosphere of constant 
thickness. Second the latitude of the most compressive meridional stress moves away from the 
equator toward the latitude where the thick zone begins. This peak in compressive stress is 
however not well-pronounced and for all practical purposes the most compressive stress can be 
said to occur in the whole thin zone. If the thin zone is very localized, stress curves become flat 
in the thick zone and approach the value predicted for a lithosphere of constant thickness (this 
is not yet apparent on Fig. [6^ because k should be much larger than 2). 

If the shell is expanding (see Fig.[T^), stresses change sign, so that Anderson's theory predicts 
normal faults striking east-west and preferably formed at the equator. Near the surface tensile 
failure is also possible because lithostatic pressure might be too small to render all stresses 
compressive. It would then lead to the appearance of joints striking east-west, because rocks are 
much weaker in extension than in compression. 

If the contracting shell is thinner at the poles, the azimuthal stress is more compressive than 
the meridional stress and they are both most compressive at the poles, so that the tectonic 
pattern consists of thrust faults striking north-south and preferably forming near the poles. If 
an expanding shell is thinner at the poles (see Fig.[T^), the azimuthal stress is more extensional 
than the meridional stress and they are both most extensional at the poles, so that the tectonic 
pattern consists of normal faults striking north-south and preferably forming near the poles. 

5.3 Despinning 

For the despinning of a shell thinner at the equator, the stresses have the following characteristics: 

• The meridional stress is everywhere extensional, whereas the azimuthal stress is extensional 
near the poles and compressive near the equator. This situation favors the development of 
normal faults near the poles and strike-slip faults near the equator. The limit between the 
faulting provinces is given by the latitudes for which the azimuthal stress vanishes. These 
boundaries are at about ±48° latitudes (independent of i^) when the elastic thickness is 
constant, and move toward the poles if the elastic thickness decreases at the equator. 

• The meridional stress is more extensional - or less compressive - than the azimuthal stress. 
This situation either favors normal faults striking east-west, or strike-slip faults striking 
at about 60° from the north (depending on whether the tangential stresses have the same 
sign or not). 

• The meridional stress is most extensional at the poles if the elastic thickness is constant. 
The maximum moves toward the equator if the elastic thickness variation is large enough 
(if a is limited to degree two, the maximum moves away from the pole if a2 < —0.44). If 
the equator-to-pole thickness ratio is small (threshold w 1/4) and the thin zone is not too 
large (threshold fc « 1 if r = 1/10), the most extensional meridional stress can be at the 
equator but it is a rather extreme situation (see Fig. EJd). 

The despinning fault pattern predicted with Anderson's theory is thus not substantially altered by 
the variation of the elastic thickness. There is an equatorial province of strike-slip faults striking 
at about 60° from the north, plus two polar provinces of normal faults striking east-west. The 
boundaries between the faulting provinces move toward the poles as the elastic thickness becomes 
thinner at the equator. As in the case of expansion, tensile failure may occur near the surface. 
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leading to the production of east-west joints. In the rather extreme case where the meridional 
stress is most extensional at the equator, these joints would first form at the equator. 

If the shell is spinning up, stresses are the same except for a sign change, so that the tectonic 
pattern consists of an equatorial province of strike-slip faults striking at about 30° from the north 
and polar provinces of thrust faults striking east-west. 

If the shell is thinner at the poles (see Fig. [TJs), the tectonic pattern is similar to the one 
for a constant elastic thickness, except that the boundaries between the tectonic provinces shift 
toward the equator as the shell becomes thinner at the poles. 



5.4 Contraction plus despinning 

Let us now consider tectonics due to simultaneous contraction and despinning. I initially assume 
that the lithosphere is thinner at the equator and I use the parameterization p3p to sta rt with 



Resul ts for the case of polar thinning are given at the end. Fig. [7] (inspired by Fig. 5 o uMelosh 



1973) represents the latitudes of the boundaries between tectonic provinces when a varying 



amount of contraction (or expansion) is added to despinning (Poisson's ratio is chosen to be 0.25 
but only weakly affects the position of the boundaries) . The proportion between contraction and 
despinning is parameterized by the contraction/despinning ratio x- 

X = -— , (57) 

which is zero if there is only despinning, positive if there is additional contraction, negative if 
there is additional expansion. Fig. |S] illustrates some possible tectonic patterns. The predicted 
tectonic pattern has the following features: 



If there is no contraction or expansion, strike-slip faults are predicted near the equator and 
east-west normal faults near the poles. The boundaries of tectonic provinces move by a 
few degrees with respect to the case of constant elastic thickness. 

If the despinning planet also contracts, the strike-slip fault province extends toward the 
poles. Beyond a first threshold of contraction, an area of north-south thrust faults appears 
near the equator; this area gets larger if contraction increases, whereas the strike-slip fault 
province is split in two smaller parts which are displaced toward the poles. Beyond a second 
threshold of contraction, the provinces of normal and strike-slip faults vanish. Thrust faults 
are then predicted over the whole planet. The two thresholds have a moderate dependence 
on the equator-to-pole thickness ratio and a weak dependence on Poisson's ratio. 

• When contraction exceeds the first threshold, thrust faults striking north-south start to 
form. Since thrust faults for pure contraction strike east-west, the orientation of the faults 
must change for a large enough contraction. If a is limited to degree two, Eqs. ([28| - ((29)) 
show that the orientation of thrust faults switches everywhere from north-south to east-west 
if W2 = a2Wo, that is if x = ^ l/a2- 

• If the despinning planet also expands, the strike-slip faulting province becomes smaller and 
vanishes for a large enough expansion. East- west normal faults are then predicted all over 
the planet. If the thickness variation is small, normal faults preferably form in polar areas. 
Beyond some expansion threshold, normal faults preferably form in the equatorial region, 
except if the thickness is constant. 

If a is not limited to degree two, the orientation of thrust faults may depend on the latitude. 
With the assumptions 52, a4 <C 1 and other cti = (i > 4), Eq. ([5^ gives the sign of the stress 



20 




Figure 7: Faulting style and orientation for simultaneous despinning and contraction (or expansion) when the 
lithosphere is thinner at the equator. The thickness variation is parameterized by the equator-to-pole thickness 
ratio r as on Fig. [2] The x-axis is the contraction/despinning ratio x defined by Eq. JSTj. Poisson's ratio is equal 
to 0.25. (a) Boundaries of the areas with a given fault style. For given 0.2 and Xi the intersections of the curves 
with a vertical line of abscissa x gi^^ the latitudes of the boundaries between tectonic provinces, (b) Contraction 
threshold beyond which the orientation of thrust faults changes from north-south to east-west. If a is limited to 
degree two as is the case here, the threshold is independent of latitude and given by x = — 1/<52. The threshold for 
r = 1 (02 = 0) is at infinity because the orientation of thrust faults never switches from north-south to east-west 
when the elastic thickness is constant. The boundaries of the areas with a given fault style are shown as in (a). 



difference at the threshold W2 = a2WQ. If < 0, the change from north-south to east-west thrust 
faults first occurs in the polar region. As contraction increases, the frontier between provinces 
of thrust faults with different orientations shifts from the pole to the equator. Conversely, the 
change from north-south to east-west thrust faults first occurs in the equatorial region if 04 > 
and the frontier moves from the equator to the pole as contraction increases. These two cases 
are illustrated by Figs. [Hh-jb where the parameterization ([55]) is used: negative values of 04 
correspond to negative values of k (enlargement of thin zone), whereas positive values of 
correspond to positive values of k (reduction of thin zone). For a given k, one can select a 
preferred latitude for the boundary between these provinces and then read on the x-axis the 
required contraction/despinning ratio %. Using the same symbols as in Eq. (j30l) . I can relate the 
contraction to the rotation variation: 

6R=~xh r \ ■ 58 

The above equation results from the substitution of Eqs. (IMl) and ([5^ into Eq. ([57]) . The 
important conclusion is that a particular combination of contraction and despinning leads to 
tectonic provinces of thrust faults differing in orientation. 

Finally I consider simultaneous contraction and despinning on a lithosphere thinner at the 
poles. In contrast with the case of lithospheric thinning at the equator, the tectonic pattern when 
contraction (resp. expansion) is dominant consists of north-south thrust faults (resp. normal 
faults) preferably formed near the poles. Hence thrust faults do not change in orientation as 
the amount of contraction increases but the area where they preferably form is displaced from 
the equator to the poles. The change in orientation occurs in the expansion regime: normal 
faults switch from striking east-west to north-south when W2 = a2'WQ (sec Eqs. with 
0:2 > 0). Near the threshold, provinces of normal faults with different orientations (north-south 
or east- west) can coexist. 
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Figure 8: Tectonic patterns predicted by Anderson's theory for simultaneous despinning and contrac- 
tion/expansion when the lithosphere is thinner at the equator. This is the pictorial interpretation of Fig. [7] 
The y-axis is the latitude. The third pattern from left is for despinning only; expansion and contraction are 
respectively added to the left and to the right of it. Normal faults, strike-slip faults and thrust faults are respec- 
tively represented with dashed, dotted and continuous lines. Faulting preferably occurs where the lines are the 
thickest. The figure does not show all possible transition states. The second pattern from the left is only realized 
if the equator-to-pole thickness ratio is close to one. The transition pattern EW-NS is shown for the case of an 
enlarged thin zone (A; < 0). 
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Figure 9: Prediction of thrust fault orientation for simultaneous despinning and contraction when the extension 
of the thin zone varies as on Fig.[3](r = 1/2; supplementary curves are drawn for k = ±1). These figures generalize 
Fig. [7]3: the vertical line in both graphs corresponds to the vertical line r = 1/2 in Fig. [TJj. Poisson's ratio is 
equal to 0.25. (a) Thin zone enlargement: latitude of the transition between the polar province of east- west faults 
and the equatorial province of north-south faults, (b) Thin zone reduction: latitude of the transition between the 
polar province of north-south faults and the equatorial province of east-west faults. 
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6 Applications 



6.1 lapetus' ridge 

Cassini images revealed in December 2004 an extraordinary feature on Saturn's sat ellite lapetus 



a hig h ridge running along the equator spanning more than half of its circumference [Porco et al 
I2OO5I I. High-resolution imaging is not available on the whole surface, but it seems that impact 



basins are present where the ridge has not been detected Denk et al. . 2008| . It is thus reasonable 



to assume that the ridge originally spanned the whole circ umference. In well - imaged areas, it h as 



an average width of 60 — 70 km and heights up to 18km Giese et al. . 2008 : Denk et al. . 200^ 



Several theories have been proposed to explain the origin of the ridge, including deposition 
of ring remnants, con vection, despinning and compaction. The morphologic analysis of the ridge 
bv I Giese et al. |2008l | ex cluded an exogenous origin due to the accr etion of a now- vanished ring 



2006j . Convection [Czechowski and Leliwa-Kovvstvnsk^ . |2008| cannot create a sufficiently 
narrow topographic rise, ev en if a width of 100 — 200 km is generously attributed to t he ridge 



Roberts and Mmmol . [2009} . Despinning [Porco et all I2OO5I: \Melosh and NimnH bOOOj cannot 



account for equatorial faults striking eas t- west, even if the lithosphe ric thickness is variable (see 
Section [531 dikes are discussed below). \Castillo-Roaez et~al ] [2007j suggested that the ridge is 
due to surface reduction caused either by despinning or by compaction, without explaining how 
this surface reduction could be localized along the equator. 

The tectonic patterns studied in this paper provide new possibilities. In particular, the 
contraction of a lithosphere thinner at the equator leads to a tectonic pattern of compressional 
faults striking east-west and preferably formed at the equator. The two assumptions underlying 
this mechanism are not arbitrary. First, the lithosphere must be thinner at the eq uator because 
of the latitudinal variation in solar insolation, as shown in Figs. 13 and 14 of \Howett et al. 
2010j . Second, a collapse of poros ity during the first 15Myr of lapetus occurs in the interior 



model of Castillo-Roaez et al. I [2007j . leading to a global contraction of about 40 km during which 
the average lithospheric thickness is less than 20 km. Though it is not known whether faulting 
occurred at the ridge, other tectonic processes such as buckling or folding will result in the same 
orientation and location since the meridional stress is more compressive than the azimuthal stress 
and the maximum compression is at the equator. 

How does despinning affect the contraction hypothesis? lapetus is at present locked in a 1:1 
resonance with Saturn: its rotation period and orbital period are both equal to 79 days. This 
synchronicity is generic of all large satellites which initially rotated much faster. Moreover lapetus 
has the peculiarity of being very flattened (equatorial and polar radii differ by 35 km), though it 
should be nearly spherical if it were hydrostatic. Its present sha pe is well fitted by a h omogeneous 
hydrostatic ellipsoid with a rotation period of about 16 hours Castillo-Roaez et 'all . [2007] . The 
most obvious explanation is that lapetus had an initial rotation period of 16 hours or less and 
that its shape froze while despinning. The lithosphere must have been several hundred kilometer 
thick when the rotation period reached 16 ho urs, otherwise the body w ould not have retained 



its flattened shape. In the evolution model oi \Castillo-Roaez et ali 2007|, despinning is a rapid 



process occurring 800 Myr after formation, at which time the lithosphere has a thickness of about 
240 km. In this scenario, contraction happens long before despinning, when the lithosphere is thin 
enough to be easily deformed. In such a case the east-west faulting pattern due to contraction 
is not directly affected by despinning. 

If contraction occurred before despinning, a question to be addressed is the absence of a 
superimposed tectonic pattern due to despinning. This problem is not unique to lapetus. Though 
all large satellites underwent despinning, none exhibits unambiguous evidence of it in its global 
tectonic grid. As for lapetus, one way out is to postulate an initial period not much shorter than 
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16 hours. If the hthosphere is thick enough, the flattening change is small and the associated 
stresses not large enough to cause faulting. 

How thick should the lith osphere be in order t o resist despinning stresses? The simplest failure 
criterion is that of Coulomb Jaeaer et which consists of a linear relation between the 



extreme principal stresses. For thrust and normal faults, this criterion can be rewritten as a 
constraint between the maximum differential stress tx/ = \(7i — (T3I (cti and (T3 are the most and 
least compressive stresses, respectively) and the vertical stress Cz — pgz: 

Of = a + b(7z , (59) 

where a is the yield strength of ice at zero pressure and CTz is positive. If < 12 MPa, 
the constants (a, &) for ice at a temperature of about lOO K are equal to (3.4 MPa, 1.86) and 
(1.2 MPa, 0.65) for compression and extension, respectively Beeman et al . 1988| . Failure starts 



either at the poles or at the equator since the differential stress is maximum at these locations. 
At the poles, normal faults start to form when the meridional stress exceeds the yield strength 
of ice for extension. At the equator, thrust faults form instead of strike-slip faults when the 
shell is thick (see Appendix 18. 4p : faulting begins when the azimuthal stress exceeds (in absolute 
value) the yield strength of ice for compression. As discussed in Section [51 a thick lithosphere can 
usually be approximated by a shell of constant thickness so that Eqs. (I108p - (|109p become valid, 
lapetus' mean radius is 735.6 k m, its mean density is 1083 kg/m^ and thus its surface gravity is 
0.22 m/s^ Thomas et al . 2007l |. Young's modulus is 11.7 GPa and Poisson's ratio is 0.325 (see 



Appendix 18. 4p . but I compute Love numbers with v = 0.5 (see Fig. (TU]). If the initial period is 
16 hours, normal faults form at the poles if the thickness is smaller than 400 km whereas thrust 
faults form at the equator if the thickness is smaller than 220 km. These thresholds are respec- 
tively reduced to 325 km and 150 km if faulting is initiated at 5 km depth. There is also a lot of 



uncertainty in the experimental value of the yield strength of ice at zero pressure: \Beeman et al. 
[l988l] consider alternative fits leading to yield strength ranges of [0, 2.8] MPa for extension and 
[0, 6.8] MPa for compression. The thickness thresholds associated with the upper values of these 
ranges are 270 km for normal faults and 125 km for thrust faults. Combining the effect of fault 
initiation at depth with high yield strength values, one can obtain a thickness threshold for nor- 
mal faults lo wer than 240 km, which is the lithospheric thickness at the time of despinning in 
the model of Castillo-Roaez et al\ 20071 ] (note that the initial rotation period is 7 hours in that 



model). Therefore, it is conceivable that lapetus' lithosphere resisted failure during despinning 
if the initial period was not much shorter than 16 hours. 

Let us now examine what happens if contraction and despinning are simultaneous. As shown 
in Figs. [7] and [HI the combination of contraction and despinning leads to three possible faulting 
patterns at the equator, given in order of increasing contraction: strike-slip faults, north-south 
thrust faults, or east-west thrust faults. The amount of contraction required to yield the third 
outcome is given by Eq. ([SO)) , at least if the lithospheric thickness varies as in Eq. (|33p . As 
a starting point, I set the initial rotation period to 16 hours, the equator-to-pole thickness 
ratio to 1/2 (i.e. 0.2 = —2/5) and the Love number to 5/2 (limit of vanishing lithospheric 
thickness). The contraction threshold between north-south and east-west thrust faults is then 
5R w 60 km which is la rge, though of the same or der of magnitude as the one due to porosity 



collapse in the model o nCastillo-Roaez et al\ [20071 ]. Of course, the amount of contraction due to 



porosity depends on the initial porosity profile which is unknown. Other values of the parameters 
strongly affect the result. On the one hand, the threshold increases with the square of the initial 
rotation rate. On the other hand the threshold decreases as the lithospheric thickness increases, 
because the flattening variation is then smaller. The two-layer incompressible model discussed in 
Appendix 18.41 vields 5R « 13 km and 6R ~ 0.9 km if the lithospheric thickness is equal to 20 km 
and 240 km, respectively. Furthermore the threshold not only decreases with the equator-to-pole 
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Figure 10: Secular degree- two tidal Love numbers for Mercury and lapetus as a function of the lithospheric 
thickness: (a) (b) ll^ /hl^ . The model is described in Appendix 18.41 For Mercury, the lithospheric thickness is 
limited by the size of the core. On an incompressible planet, despinning leads to north-south thrust faults in the 
equatorial region when l"^ /h"^ < 1/4 (see Eq. llllOI I with u = 0.5); this threshold is indicated by the horizontal 
line in graph (b). 



thickness ratio (see Fig. [7}d), but is also affected by the precise form of the thickness variation: 
Fig. |nt> shows that a reduction of the thin zone can decrease the value of the threshold at the 
equator by a factor of two. Therefore, simultaneous contraction and despinning may lead, or not, 
to east-west thrust faults. For a given amount of contraction, the outcome sensitively depends 
on the initial rotation rate and on the lithospheric thickness (average thickness, equator-to-pole 
thickness ratio, thin zone size). 

Besides the contraction hypothesis, an episode of expansion in lapetus' evolution due to 
differentiation Sauvres and Croft . 1986{ would lead to equatorial normal faults striking east- 
west and forming first at the equator. lapetus' ridge does not look like normal faulting, but the 
sa me extensional stress field i s compatible with the intrusive dike origin for the ridge proposed 
by \Melosh and Nimmo 2009| . In that scenario, the simultaneous occurrence of despinning is 



generally not a problem because it is only required that the meridional stress be most extensional 
at the equator. This is true for moderate amounts of expansion, the value of the threshold 
depending on the thickness profile. The expansion threshold is for example 7.5 times smaller in 
magnitude than the contraction threshold discussed above, assuming an equator-to-pole thickness 
ratio of 1/2 (it coincides in this particular case with the expansion threshold for normal faults at 
the equator). Note that despinning can also be favorable to the formation of an east- west dike 
at the equator, but the conditions are rather stringent: the equator-to-pole thickness ratio must 
be smaller than 1/4 and the thin zone not too large (see Fig. [SJj). 



6.2 Mercury's scarps 

Lobate scarps on Mercury are linear or arcuate escarpments varying in length from about 



20 km to 600 km and in height from a few hundred meters to about 1.5 km [Strom et al. ■ 11975 
Watters et al. , 2009| . They are the most common tectonic features on Mercury and are thought 



to be the result of thrust faulting. The predominant north-so uth orientation of lobate scarps was 
first attributed to illumination bias but \ Watters et all |2004| showed that this preferred orienta- 



tion is real. Th e combination of planetary contraction and despinn i ng is a possible explanation 
for th is pattern Melosh and Dzurisin . 119781: \Pechmann and Melosm . Il979t iDombard and HaucM . 
2008| . Another interesting observation is that lo bate scarps in the sou thern polar region are as 
often east-west than north-south (see Fig. 2.23 in Watters and Nimm^ 2009.] ). East- west lobate 
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scarps are much rarer in the northern polar region, but the total length of scarps of any orien- 
tation is also much smaller than in the south. Let us suppose that lobate scarps predominantly 
strike north -south from 50° S to 50° N latitude and east-west in the polar regions, as done by 
2008l |. Fig. shows that simultaneous contraction and despinning may lead to this pat- 



Kino 



tern if the thin zone is rather extended. As mentioned in the introduction, th e lithosphere could 



be thinner by a factor of two at the equator in comparison with the poles [McKinnon . 11981 



Melosh and McKinnon . 1988| . In that case, the required contraction/despinning ratio should 
be close to two, assuming that the shape parameter k is between —1 and —2. Substituting 
X = 2 in Eq. (j58p . I obtain SR ~ 20km if h"^ = 5/2 and the initial rotation period is 20 hours 



{g = 3.7 m/s^, R — 2440 km). The contraction threshold is somewhat reduced by more realistic 
values of hj . The thickness of the lithosphere at the time of the formation of lobate scarps 
is difficult to pinp oint: arguments have been made for values of 25 — 30 km and 100 — 200 km 
(see a summary in I Waiters and Nimmo [2003). In any case, Fig. [TUk shows that the secular 
Love number is not very sensitive to the lithospheric thickness. The contraction threshold 
is only reduced to 15 km and 12 km if the lithospheric thickness is equal to 30 km and 100 km, 
respectively (see Fig. [TUl and Appendix | 8.4I for the value of hT^). The radius contraction at the 



origin of lobate scarps likely being in the range 1—2 km Strom et al. . 1975 : Watters et al. . 2009j , 



it is far too small to make possible the existence of east-west thrust faults. 

One kilometer of contraction, combined with despinning, is also too small to generate north- 
south thrust faults all over the planet, though two kilometers are enough if the yield strength of 
rock is neglected (use Eq. ((55)) with x determined from Fig. [7^). Yet another problem is that 
lobate scarps are more abunda nt in the south than near th e equator, contrary to the prediction 
of contraction plus despinning. Dombard and Hauc^ 2008| thus suggested that a larger amount 
of contraction occurred during despinning in the early history of Mercury (their contraction 
threshold of 3 — 5.5 km is however significantly reduced if the factor is taken into account). 
The Late Heavy Bombardment then erased the faults and new scarps were produced along the 
old fault lines by a later global contraction event of smaller magnitude, when the planet had 
despun. Such a scenario widens the choice of the relative weights of despinning and contraction. 
On the one hand despinning needs not be finished when the early tectonic grid forms. On 
th e other hand early contract ion can be rather large in thermal evolution models: 2.5 — 7.5 km 



in \Dombard and HauciA |2008j |. with an upper bound of 17 km if the core completely solidifies 
[Solomorl '. 1976| (me asurements of Mercu ry's librations however indicate that the core is at least 
partially molten, see \Maraot et all 2007 |V In conclusion, it is possible to imagine a sequence 
of events in which various events of contraction and despinning lead to thrust faults that strike 
north-south near the equator and cast-west near the poles. However the observation that the 
total scarp length increases from north to south remains to be explained. East-west thrust 
faults could instead be the result of the reactivation of early normal faults due to despinning 
Watters and Nimmd l2009l | . The various orientations of loba te scarps could also be due to a 



combination of despinning, contraction and true polar wander Matsuvama and Nimmo . 2009| . 



7 Conclusions 

The main result of the paper can be very simply stated: the contraction or expansion of a 
planet with lithospheric thinning at the equator results in tectonic features striking east- west, 
preferably formed in the equatorial region. If the lithosphere is thinner at the poles, contraction 
or expansion generates tectonic features striking north-south. In general, despinning alone cannot 
produce east-west features near the equator. One exception could be the appearance of east- 
west joints or dikes at the equator if the thickness variation is large and strongly peaked at the 
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equator. 

A spectacular and rather unique illustration of an east-west structure is the equatorial ridge 
on lapetus, which can be explained by the theory developed in this paper. The mountainous 
structure of the ridge is suggestive of a compressional event but an extensional process is not 
excluded. Since despinning inhibits the formation of east-west compressional tectonics, the ridge 
probably formed long before lapetus had despun, but this is not required if the contraction is 
large enough. 

Contraction and despinning stresses can be computed either numerically or with a semi- 
analytical method, which yields a perturbation expansion in the parameters describing the thick- 
ness variation. First order formulas for stresses are sufficient if the thickness varies by a factor 
of two. In the simplest case (inverse thickness of degree two) , contraction stresses are well ap- 
proximated by the elementary formulas while Eqs. (jl28p - (|129p are valid at first order 
for an arbitrary thickness profile. 

The combination of radius change and despinning generates new tectonic patterns when the 
lithospheric thickness is variable. On the one hand, contraction plus despinning makes possible 
the coexistence of tectonic provinces of thrust faults differing in orientation if the lithosphere 
is thinner at the equator. Lobate scarps on Mercury roughly follow such a pattern but the 
small amount of contraction associated with the scarps makes it necessary to resort to rather 
complicated models involving fault reactivation. On the other hand, expansion plus despinning 
may lead to the coexistence of tectonic provinces of normal faults differing in orientation if the 
lithosphere is thinner at the poles. In both cases the magnitude of the contraction/despinning 
ratio must be large otherwise the predicted patterns are similar to those valid for a lithosphere 
of constant thickness. 

Global contraction or expansion occur in many models of interior evolution and are often 
assumed to be the cause of planetary tectonics. Since lithospheric thinning due to the latitudinal 
variation in solar insolation must be a generic phenomenon, it is surprising that an east- west 
orientation of tectonic features at the planetary scale is so rare in the solar system. One reason 
could be that the thinning is too small to have a significant effect. Also, tidal heating might 
counteract the effect of the variation in solar insolation. Another explanation is the simultaneous 
occurrence of despinning. Finally, it is possible that global contraction or expansion events 
generally happen very early in the history of the planet so that either faults do not appear 
because the lithosphere is not yet formed, or faults do form but the associated tectonic pattern 
is subsequently erased. Further inferences about the origin of tectonic features depend on more 
complete tectonic mapping, which is underway for Mercury with Messenger data and for Saturn's 
icy satellites with Cassini data. 
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8 Appendices 



8.1 Differential operators on the sphere 

The operators Oi are linear differential operators of the second degree on the sphere: 



O3 



80^ 



cot ( 



1, 

\dddLp dtp 



d_ 

de 



1, 



(60) 
(61) 
(62) 



where 9 is the colatitude and Lp is the longitude. They give zero when applied on spherical surface 
harmonics of degree one. 



The scalar differential operator A' is defined by \Beuthd . l2008j : 



A' 



= O1 + O2 
= A + 2, 



(63) 
(64) 



where A is the spherical Laplacian (called surface Laplacian in BeuthA 2008llV Spherical surface 
harmonics of degree i are eigenfunctions of A' with eigenvalues Blakelii . 1995 . pp. 121-122] 



S'^ = -i{e + 1) + 2 . (65) 
The scalar differential operators C and A appearing in the membrane equation ^ are defined 



by 



C(a;&) = A' (a A' 6), 

A{a;b) - {Oia){02b) + {02a){Oib)-2{03a){03b). 



(66) 
(67) 



Since C{a;b), A{a]b) and A'b give zero if & is a spherical surface harmonic of degree one, the 
membrane equation ^ does not constrain the degree-one components of the stress function F 
and of the transverse displacement w (translation invariance). 

C and A are linear in a and in b. If a is constant A{a ; b) = aA'b whereas A{a ; b) = bA'a if b 
is constant, so that the following identities hold: 



C{l;b)-{l + v)A{l;b) 
C{a;l)- {l + v)Aia;l) 



(A' - 1 - i^) A'5 . 
(l-jy) A'a. 



(68) 
(69) 



8.2 Operators Oi on Legendre polynomials 



Legen dre polynomials Pi (x) are defined in many books [Whittaker and Watsom . Il935t \BlakeM . 
Il995l |. I only give those that appear in the text, in the form Pi{cos6): 



P2 
Pi 



= 1 



(3 cos 261 + 1) , 



64 
1 

512 



(35 cos 46* + 20 cos 26* + 9) 



(231 cos 661 + 126 cos 461 + 105 cos 26* + 50) 



(70) 
(71) 

(72) 

(73) 
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The action of the operators Oi defined by Eqs. (|60p - (p^ on Legendre polynomials Pe{cos9) 

is 



OiPe 
O2P1 
O3P1 



(f - Pi + 

(f-^)P,-E£, 

0, 



(74) 
(75) 
(76) 



where 



r=l 



= (2^-3)P,_2 



1) Pl^2r , 

- (2£ - 7)P,_4 



(77) 



with p — £/2 for even £ and p = {£ — l)/2 for odd £. The action of O2 is obtained from 
an ex pression for the first derivative of a Legendre polynomial found in I Whittaker and Watson 
I935I p. 330]. The action of d is then obtained from Eqs. §^ and ([SS]) . 



The values of OiP^ and O2P1 are equal at the poles: 

1 



O^Pe\e=o = 



O^PAe=^ = 



2 ^ 
(-1)^ 



(« = 1,2), 
S'i (i = l,2) 



(78) 
(79) 



where S'^ is defined by Eq. (I55|) . 

Since the coefficients multiplying the Legendre polynomials in Eqs. (|74p - ([75)) are independent 
of 9, the operators Oi can be represented as matrices acting on vectors of Legendre coefficients. 
The action of Oi on a vector of Legendre coefficients of degrees 0,2,4 and 6 is represented by 
the matrix 

/I 1 1 1\ 



(6) 



(80) 



0-3 5 5 
-15 9 

V -35 / 

The action of O2 on a vector of Legendre coefficients of degrees 0,2,4 and 6 is represented by 
the matrix 

/I -1 -1 -1\ 

(6) _ 0-1-5 -5 
0-3-9 

\ -5 / 



O 



(81) 



8.3 Constant thickness: thin shell 

When the shell thickness h is constant, the membrane equation ([2]) becomes 

(A' - 1 - !/) A'F =^A'w, 
R 



(82) 



where I used the identity (j68p . Eq. (j82p is diagonalized by expanding F and w in spherical 
surface harmonics whose coefficients (£ 7^ 1) satisfy (see Eq. (1551) ) 



Fin 



Eh 



1 



R -£{£+!) + 1-V 



wen 



(83) 
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The index m refers to the order which is equal to zero under the assumption of axial symmetry. 
The component of degree one remains indeterminate. The stress function generates the stresses 
according to Eqs. (H])-®. The action of Oi on Legendre polynomials of degree 0, 2, 4, 6 is given 
by Eqs. (|80|-(|8T1). 

For a contracting shell, w — wqPo and the stress function is constant (Pq — ^)'- 



The non-zero stresses are 



pC ^ Eh wq 

l-iy R °' 

E Wo 



c ^ 



\-v R 

For a despinning shell, w — W2P2 and the stress function reads 

Eh W2 „ 



5 + 1/ R 



The non-zero stresses are 



n 



n 



E W2 



b + v R 



{-P2 1) , 



(84) 
(85) 

(86) 
(87) 



h + v R 

The coefficient W2 can be related to the flattening variation 5f {Sf < for despinning) by 

2 



■W2 



RSf. 



(89) 



After substitution of Eqs. ([7T|) and despinning stresses Venina-Meinesz . 1947 : MelosfA . 



mm read 



"98 

n 



E 5f 

4 (3 cos 261 -I- 5) 



The flattening is related to the angular rate by 



5 + 6 

E Sf 

4 (9cos26i- 1) . 



(90) 
(91) 

(92) 



where q is the surface gravity an d /i^ is the degree-two tidal Love number for radial displacement 



Matsuvama and Nimmd . [2008 1. 



For a deformation of degree four, w — W4P4 and the stress function reads 

Eh 
19 + 1/ "R 



The non-zero stresses are 



^ ^(-3P4-5P.-1) 



19 + i/ P 
E W4 
19 + 1/ P" 



(-I5P4 + 5P2 + I) . 



(93) 

(94) 
(95) 
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8.4 Constant thickness: thick shell 



In a tidal or centrifugal potential, the radial deformation of the surface of a spherically stratified 
planet is parameterized by the Love number /il", as in Eq. (1921) . Other Love numbers noted I2 and 
k'^ parameterize the tangential deformation of t he surface and the gravity perturbation at the 



surface induced by a tidal or centrifugal potential \Lambeck . 1980| |. Love numbers are complicated 
functions of the internal planetary structure. A simple formula exists for an incompressible body 
of uniform density and shear modulus (see below), but this model is not appropriate for long-term 
processes such as despinning. Stresses indeed relax by viscous flow or creep everywhere within 
the planet except in the lithosphere, which is by definition rigid on geologic time scales. Internal 
layers in which stresses have relaxed are modeled as flui d, inertial forces are neglected (static 
limit) and the resulting Love numbers are called secular Matsuvama and Nimmo . 2008| . The 
simplest interior model including a lithosphere is a two-layer spherically stratified planet made o f 
an incompressible fluid core surrounded by an incompressible elastic lithosphere [MelosM . Il977 1 . 
A more realistic density profile is obtained by dividing the planet into several homogeneous 
layers, each of constant thickness. Tida l Love numbers for such models are readily computed 
with the propagator matrix technique [Sabadini and Vermeersera . |2004| . The assumption of 
incompressibility is not realistic but satisfactory for our purpose, since the main uncertainty 
regarding secular Love numbers lies in the unknown thickness of the lithosphere. 

In the limit of vanishing core radius, the tw o- layer mode l gives the degree- two tidal Love 
numbers for a homogeneous incompressible body \Love . 1944 . pp. 257-259]: 



hi = 



1 



19 



2 pgR 



10 ■ 

- h'^ 



(96) 

(97) 
(98) 



R is the radius, p is the dens ity, g is the surface gravity and /x is t he shear modulu s. Note that 
Munk and MacDonald 1960| give an incorrect formula for while Lambeck 1980| gets it right. 
Love numbers for a hydrostatic homogeneous planet are obtained by setting /i = 0: = 5/2 
and fc|" = 3/2. The Love number I2 does not exist for a hydrostatic body since displacements in 
a fluid are indeterminate in the static limit, apart from the radial displacement at the surface. 

In the limit of vanishing elastic thickness, the two-layer incompressible model gives the fol- 
lowing degree-two tidal Love numbers: 



5 

2 ' 

11^2 

- 

5 2 ■ 



(99) 
(100) 
(101) 



As expected, and fc^" coincide with the values for a homogeneous hydrostatic body. More 
generally, a lithosphere of vanishing thickness neither affects the shape of the body nor its 
gravitational field, whatever the internal density distribution. However I2 differs from the value 
for a homogeneous incompressible body since the presence of a membrane modifies tangential 
displacements at the surface. 

I now show that the value of given by Eq. (llOOp agrees with the membrane limit of 
thin shell theory. Toroidal displacement is absent when the shell thickness is constant (axial 
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symmetry is actually sufficient), assuming that there is no toroidal external potenti al. The 



tange ntial displacement v can then be expressed as the surface gradient of a potential S \Beuthe , 
I2nn8l| : 

w^VS. (102) 

The displacement potential S is related to the stress function by Eq. (71) of Beuthe 2008| . which 
becomes in the membrane limit 

lis ^ Ra{l-v)MF ~2w . (103) 

The substitution in this equation of the stress function yields the degree-^ component of S 
in terms of wf. 

^^-^(TTTFTT^-- (104) 

If V2 = — ri^i?^/3 is the degree-two component of the external potential at the surface, S2 and 
W2 are related to degree-two tidal Love numbers by 

W2 = —V2, (105) 
.9 
;T 

^2 = —V2. (106) 
9 

In the membrane limit, thin shell theory thus predicts that 

iT^lUih^, (107) 

This relation holds for any density distribution within the planet, as long as it is spherically 
symmetric. If the lithosphere is incompressible (v = 1/2), the ratio I2 /h'^ is equal to 3/11, that 
is Eq. (jlOOp . If the lith osphere is made of ice {v w 1/3), the same ratio is equal to 1/4 (this 
result has been noted by ' Wahr et al) '20095 regarding tidal stresses). 

The insertion of Eq. ((89)) and Eqs. (I105p - (ll06p into stress-strain relations yields the following 
expressions for the tangential stresses at the surface of a spherically stratified planet: 

^Fe = -j^'^5f[{l + v)P2+^^{l-v-2{2 + v)P2)^, (108) 

<^ = - \ ((1 + ^2 + ;^ (-1 + - 2(1 + 2i.)P2)) , (109) 

where I used the fact that the rr component of the stress vanishes at the surface. In the 
membrane limit of thin shell theory. Love numbers satisfy Eq. (|107p so that Eqs. (|108l) - (|109p 
become equivalent to the stresses Mo reover the subs titution of Eq. (1^ into Eqs. (llOSp - 

(|109p yields the same stresses as those found hy lMelosf} 1977| for a homogeneous incompressible 
body (see his Eq. (20)). 

Once Love numbers are known for some interior model, it becomes possible to compute surface 
stresses and thus predict despinning tectonic patterns using Anderson's theory of faulting (see 
Section [5]). This is a conv enient way to compute thick shell effects on faulting style as was 
first done by \Melosh 1977 1, since several codes for the computation of Love numbers are now 
available. The boundaries of tectonic provinces depend only on two numbers: the ratio Zj/ft.^ 
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and Poisson's ratio v. For example, despinning induces an equatorial zone of north-south thrust 
faults if the meridional stress is compressional at the equator, that is if 

This inequality is satisfied neither by thin shell theory (see Eq. (jl07l) ) nor by the homogeneous 
incompressible model (see Eq. (j97p ). but is verified for the incompressible models shown on 
Fig. [TUb if the shell thickness is between 8.4% and 62% of the radius. Though and I2 are 
very sensitive to the internal structure (see Fig. ITUk). their ratio varies much less. That is why 
the despinning tectonic pattern is not much affected by the shell thickness. 

For lapetus, I adopt an incompressible t wo-layer mode l with a radius of 735.6 km and a 
uniform density of 1083 kg/m^ in both layers Thomas et aZ.I . [2007 |. The shear modulus /x of ice 



at 257 K is 3.52 GPa so th at its value at 80 K (a. pproximately the mean temperature of lapetus' 
surface) is about 4.4 GPa Gammon et aZ.I . [l983{ . Poisson's ratio v for ice is 0.325 (but v = 0.5 



in the model) and Young's modulus E ~ 2(1 + v)^ is 11.7 GPa. The model yields = 0.53 and 
= 0.038 when the lithospheric thickness is equal to 20 km and 240 km, respectively (for other 
values, see Fig. [TUl). 

For Mercury, I adopt an incompressible three-layer model with a fluid core, fluid mantle 
and elastic lithosphere. P arameters are chosen to be in the middle range of admissible val- 



ues [Rivoldini et all . I2009l | : core radius of 2000 km, density of mantle and lithosphere equal to 
3100 kg/m'^ (thus core density equal to 7325 kg/m^), seismic wave velocities in the lithosphere 
given by Vp = 7900 m/s and Vg ~ 4550 m/s. The corresponding elastic parameters for the litho- 
sphere are: ^ = 64 GPa, E = 160 GPa, v = 0.25 (but 1/ = 0.5 in the model). The model 
yields /i^ = 1.83 and = 1.49 when the lithospheric thickness is equal to 30 km and 100 km, 
respectively (for other values, see Fig. [TU|) . 

8.5 Two properties of the strain 

The tangential displacement in the azimuthal direction vanishes if there is axial symmetry, i.e. 
if the elastic thickness and the radial deformatio n of the shell only depend on the latitude. In 



that case, the strain-displacement equations [e.g. \Beuthd . l2008l | read 

1 f dvg , 



R \ de 



(111) 



= ^{cot9ve + w), (112) 
H 

eev = 0, (113) 

where w is the radial displacement and vg the tangential displacement in the meridional direction. 

By symmetry, vg vanishes at the poles. The average of egg on [0, vr] is thus related to the 
average oi w — w/R: 

<egg>^<w>, (114) 

where 

</>=- r fde. (115) 

I" Jo 

There is equatorial symmetry if the elastic thickness and the radial deformation are symmetric 
about the equatorial plane. In that case, vg vanishes at the equator so that the average can be 
done on [0,7r/2]. For contraction plus despinning, w — iJuo + 102^2, so that 

< egg > ^ Wo + — ■ (property 1) (116) 
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Another consequence of axial symmetry is that the azimuthal strain becomes independent of 
the elastic thickness at the equator: 



ewl9=f=^le=f (117) 

This property also directly results from the variation in length of the equator, 6L — 2TTw\g=iL, 
which is uniformly distributed along the equator. For contraction plus despinning, 



£vvle=f =wo - ^ ■ (property 2) (118) 



8.6 Contraction solution at first order 

Given the despinning solution on a shell of constant thickness, the dualities of Section [5^ generate 
a first approximation of the contraction solution on a shell of variable thickness. If a depends only 
on the parameter a2 as in Eq. (j20l) . the contraction and despinning solutions can be expanded 
in 0(2 about the solution for a shell of constant thickness given in Appendix 18.31 The dualities 



(PT|l - (|24p then allow to compute the contraction solution at order aj^^ in terms of the despinning 
solution at order a2 ■ The nondimensional despinning stress function at zeroth order is obtained 
from Eq. ([86]): 

^--T^ + 0(«2). (119) 

W2 + ly 

Nondimensionalization has been carried out with Eqs. (in))- ([T^ . The duality (|21l) yields the 
contraction stress function at first order in 012: 

^^J—+a2^ + 0{al). (120) 

The duality (^5]) acting on Eqs. (|57| - (|55|) yields the contraction stresses at first order in 02 
(an alternative is to use Eqs. (fT5|) - (fro|) on F*^): 

-^^e ^ 1 2{2 + u)P2-l + iy 

— = (i_,)(5 + ,) ' (121) 

^ 1 _ 2(1 + 2t.)P2 + 1 - ^ 

- = 'a +"2 T. TTF"^ ^ • (122) 

Fig. [T5b compares Eqs. (jl21l) - (|122p with the stresses computed at high order in 0:2, that is with 
a truncation degree n — 20. The first order approximation is good for small 012 (0:2 = ^2/5) but 
is bad for large a2 (0^2 = —6/7). 

The strains at first order in 0,2 can be computed with Eqs. (jl7l) - (jl9p : 

^ = 1 + «2^(4P2-1), (123) 

1 + a2 \^ {2P2 + 1) . (124) 



Wo 



With the values P2|e=o = 1, P2\e=TT/2 = ^1/2 and < P2 >= 1/4, one can check that the stresses 
and strains (|12ip - (|124D satisfy the properties enumerated at the beginning of Section lXTl Besides 
the meridional strain as given by Eq. (|123p is equal to its average value at 6* = 45° (where 
P2 ~ 1/4), so that the meridional strain curves for varying 0:2 seem to have a unique crossing 
point (see Fig. [T^). Since this property is only valid at first order in a2, I call this point a 
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pseudo-node. Pseudo-nodes for stresses are obtained by setting to zero the term of order a2 in 
Eqs. p^ - p^ . yielding respectively 9 = 48.2° and 9 = 65.9° (j/ = 0.25) for the meridional and 
azimuthal stresses (see Fig. [S^). 

The same method works for a not limited to degree two. The generalization of the duality 
(PT|) to an a depending on an arbitrary number of parameters ai reads 

^ = 1 E"i^' 125 

where F^ denotes the solution of the membrane equation (|14p for w = w^Pl. At first order in 
ae, the contraction stress function can thus be computed from the solutions (|83p on a shell of 
constant thickness: 

where 

fe, ^ £{£+!) + (127) 

There is no contribution of ai to the contraction stress function (jl26p . but this is only valid at 
first order in ai (terms in a" vanish, but mixed terms like aia2 do not). At first order in ai, 
the contraction stresses read 



cr£, - , 1 



Wc) 1 — v ^ 



^ci, (^ + — (128) 



02PI 


Pi 




1 — V 








1 — V 



Wo 1 — v 

with the action of Oi on Pg defined by Eqs. ([74 | - ((75)) . The sum now starts sA t — 1 because 
of the multiplication by a when computing the stresses from the stress resultants, but the first 
term in the parentheses vanishes for £ = 1. At first order in a^, the contraction strains read 

# - 1-(1 + ^)^|1(0,-1)P,, (130) 



e>i 



1-(1 + '^)Et^(^2-1)P.. (131) 



Wo ^ L 

The stresses (|128p - (ll29p and the strains (|130p - (|131l) satisfy the properties enumerated at the 
beginning of Section [OI because of Eqs. ([75|) - d7^ and since Eqs. (I5n|) - (I5T]) with axial symmetry 
imply 

<{Oi-l)Pe> = 0, (132) 
i02-l)Pi\e=f = 0. (133) 



8.7 Operators C and A on spherical surface harmonics 

If the shell thickness is constant, the membrane equation reduces to Eq. (j82p which depends only 
on the operator A' and can be solved with spherical surface harmonics as in Appendix 18.31 If 
the thickness is variable, the membrane equation is not diagonal in the basis of spherical surface 
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harmonics. This basis remains however useful if the action of the operators C and A on spherical 
surface harmonics generates a linear combination of the same functions, so that the operators 
can be represented as square matrices. The action of the operators Oi - from which C and A 
are built - on Legendre polynomials is known (see Eqs. (l74l) - ((76t ) and could be used to derive 
a formula for the action of C and A on Legendre polynomials. I will however obtain a more 
compact formula for the action of A by rewriting it exclus ively in terms of the scalar operator 
A'. The operator A can be written in the equivalent form \Beuth^ . 12008 1: 



A{a-h) = (A a) (A 6) - {ViV j a){V'V^ h) + {Aa)b + a{Ab) + 2ab, (134) 

where is the covariant derivative on the sphere and summation on repeated indices is im- 
plicit (indices are raised with the inverse metric g'^^ = diag{l, sin~'^ 9)). Repeated covariant 
differentiation yields the following identities: 

A{ab) = {Aa)b + 2{V,a){V'b)+a{Ab) (135) 
AA(a6) ^ iAAa)b + 2{Aa){Ab) + a{AAb)+ AiV.,Vja){V'V^b) 

+ 4 (V,A a)(V*&) + 4 (V, a){W'Ab) + 4 (V,a)(V'6) . (136) 

In the derivation of the last identity, I used the commutation relation for the spherical Laplacian 
and the covariant derivative on the surface of the sphere, 

AV,a- V,Aa = Vjfl, (137) 

which can be de rived from the general commutation relation of covariant derivatives of a vector 
[see Eq. (Gl) in ' BeuthA l2008l |. Other useful identities are obtained by substituting a — Aa or 
& — T' A & in Eq. (I135p . With the above identities, the operator A can be rewritten as an expression 
where the only intervening operator is A': 

A{a;b) = ^[-A'A'{ab)~{A'A'a)b~a{A'A'b) 

+ 2 (A' a)(A' 6) + 2 A'((A' a)b + a (A' b)) 
-2A'{ab)-2{A' a)b-2a{A'b) + S,ab] . (138) 

The action of the operators C and A can be partially evaluated if a and b are spherical surface 
harmonics of degree m and n (noted and 6„): 

C(a„;5„) = S',,A' (ambn), (139) 
A{am;bn) = -i(A'A' + «;™„A' + A„™)(a™6„). (140) 

where 5^ is given by Eq. ((65)) and 

Kmn = 2{1-S'„^-S'J , (141) 

A,„„ - {6'„,-S'S + 2{S'„,+6'„)-8. (142) 

The action of the operators C and A can only be fully evaluated if the product a„i&„ is expanded 
in spherical surface harmonics. 

8.8 Operators C and A on Legendre polynomials 

The action of the operators C and A on spherical surface harmonics and 6„ is given by 
Eqs. (I139p - (|140p in which products b„ appear. I am now going to evaluate these products 
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with the assumption that the spherical surface harmonics are zonal or, equivalently, Legendre 
polynomials. 

The product of two Legen dre polynomials can be expand ed into Legendre polynomials with 
the Adams-Neumann formula Whittaker and Watsoru . ll935L p. 331]: 



-Pm (•^) (•^) — ^ ^ ^mnj -fm+n — 2j (•^) 
j=0 



with the coefficients of the expansion being given by 



h ■ — 

m + n + i — J 



where 



(2m -1)!! 1.3...(2to-1) 



(143) 



(144) 



(145) 



m! 771! 

with Kq — 1. These coefficients kmnj are symmetric in (771, n) and fc,„oo — ^Omo = 1- 

The action of the operators C and A on Legendre polynomials of degree 777, and ti is thus (the 
dependence on x is omitted): 



min{m,n) 
mi7i{m,n) 



' \mnj Pm-\-n — 2j ; 



mnj '4^mnj ^m-\-n — 2j i 



P„ 



where 



A' A' 

Xmnj — 0„_|_„_2j 



(146) 

(147) 

(148) 
(149) 



with (5^ given by Eq. ([55]) and coefficients (Km„, A„i„) given by Eqs. (|14ip - (|142p . If 777 = 
(constant clastic thickness), the only relevant coefficients are 



XOnO — {Ki) 
V'OnO = 5'^ ■ 



(150) 
(151) 



The Legendre polynomial of degree one does not belong to the images of the operators C and 
A: if 777, + 77, — 2j = 1 (meaning that [777 — 77] — 1), then S'^^„_2j = and A„m = 0, so that 
Xmnj = 'ipmnj — 0. The more general proof th at spher i cal su rface harmonics of degree one are 



not included in the images of C and A is given i mBeutha (2008|. Since the coefficients multiplying 
the Legendre polynomials in Eqs. (jl46p - (ll47p are independent of 6, the operators C and A can 
be represented as matrices acting on vectors of Legendre coefficients. 



8.9 Membrane matrix for n = 6 

The membrane matrix formalism described in Section|3]is illustrated here with a simple example: 
a is limited to degree two and the truncation degree 77 is equal to 6. The action of the operator C 
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on the vector F = (Jq, F2, F4, Fe)'^ can be computed with Eq. (jl46p . The matrix C'^^ appearing 
in Eq. ([52)) is given by the sum C'^^ = C, 



(6) 




where 



/ 4fcooo 






V 





~8 fc200 








16A;o20 



-8 fc222 

16A:22i 
72A:220 







324fco4o 




72 k242 

324 fc24i 

720 fc240 



1600fco60 / 






720^262 
1600^261 / 



(152) 



(153) 



The action of the operator A on the vector F is computed with Eq. (I147p . The matrix A^^^ 
appearing in Eq. (j52p is given by the sum A'^' ~ A, 



(6) 



«2Af , 



where 



A (6) _ 



2 kooo 











-4 fco20 











~18fco40 














-4^222 





-4fc200 


20 fc221 


6 ^242 





6fc220 


90 fc241 








24 fc240 







-40fco6o J 



(154) 



(155) 





(6) _ -4fc200 ^Ufc221 0^242 

24fc262 

V 24fc240 200fc261 J 

The four rows are associated with the projection of the membrane equation on harmonic degrees 
0, 2, 4, 6. The upper index on the matrices denotes that the expansion is hmited to degree 6. The 
integer coefficients are generated by the functions Xmnj and ^pmnj defined by Eqs. (I148p - (|149p . 
In particular, the integer coefficients in the matrices Cq^-* and Aq^' are given by Eqs. (jl50p - (jl51l) . 
The Adams- Neumann coefficients kmnj-, given by Eq. (|144p . are in symbohc form so that their 
origin is apparent: the mn indices indicate that the term comes from the coefficients of a and 
F of degree m and n, respectively, whereas rn+n— 2j is equal to the harmonic degree associated 
with the row of the matrix. Apart from the trivial values fc,„oo = fco?riO = 1, the values relevant 
to the example above are associated with the products P2P2, P2PA and i-2^6- 



(^220,^221,^222) = 

(^240,^241,^242) = 
(^260,^261,^262) = 



/18 2 1 
V 35' 7' 5 



5 20 2 
IT' 77' 7 
28 14 45 
65' 55' 143 



(156) 
(157) 
(158) 



The coefficient ^250 does not appear in Cj^^ and Aj^'' because it is associated with a polynomial 
of degree 8. 

The approximate membrane equation (1511) now takes the following form: 



M'^' {Fq, F2, Fi, Fef = {2wo,-4iu2,0,0f 



(159) 
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The coefficients Ff of tfie contraction solution up to a2 are generated with the recurrence 
relation starting with F(°) = {wo/{l - v), 0, 0, 0): 



^ 1 2 1 _2 4 1 + 5j/ 



Wo l-v 5 5 + ^ 35 (5 + z/)2 

F^' _ 1 _ 2 1 + 5;/ , 12 49 - 7z/ + 17t/- 



(160) 



Q2 + ^ , .,A2 "2 + — , I .A "2. (161) 

(162) 



Wo 5 + 1/ 7 (5 + z/)2 ^ 7 (5 + 1^)3(19 + 1/) 

Ff _ 6 i/-ll _2 36 (;/- ll)(-21 + 56i/ + 5z/2) 



Wo 35 (5 + i/)(19 + i/) ^ 385 (5 + z/)2(19 + z/)2 

at (163) 



^ 18 (i^-ll)(z.-29) _3 



Wo 385 (5 + I/) (19 + i/)(41 + ;/) 

Coefficients F/' have a leading term of order a2 ■ Thus the truncation at harmonic degree n 
can be understood as an expansion in 0:2 : harmonic coefficients that are neglected are at most of 
order cij^"^^: if n = 6, Ff with £ > 8 are at most of order ©(al) and are ignored. However the 
truncation at degree n is not equivalent to an expansion at order cij^^ since the lower the degree 
of the harmonic coefficient, the better it is known. More precisely, the band-diagonal structure 
of the membrane matrix has the effect that, for a given truncation degree n, the expansion of 
a term Ff can be done up to order cij ^^^^^^ included: if n = 6, {Fq , F^ , F^ , Ff) can be 
expanded up to order (0:2 ) respectively. 

The coefficients Ff of the despinning solution are related to the coefficients of the contraction 
solution by Eq. (pT|) : 



F,^ _ 1 (F§ 



W2 02 V Wo I — ly 



(164) 



F^ 1 pC 

= (i = 2,4,6). (165) 

W2 0L2 Wo 

Coefficients Ff thus have a leading term of order 03^^ ^, except Ff which begins at order 
0:2- Furthermore, for a given truncation degree n, coefficients Ff can be expanded up to order 

-n-i/2 
"2 

Figs. [TT] and [T^ show the effect of the truncation degree n on the coefficients of the contraction 
solution. The effects on the despinning solution are similar since they are related by Eqs. (|164p - 
(|165p . Since the stresses are the quantities of physical interest, the figures do not show the 
Legendre coefficients of the stress function F but rather those of the sum of the nondimensional 
stress resultants: 

m = (Ngg + N^^)^ ^S'f Fi. (166) 

Truncation effects can be large close to a2 = 2, but I only consider negative values of a2 relevant 
to equatorial thinning. Fig. [TT] shows the effect of the truncation degree and of the expansion in 
a2 on the first four Legendre coefficients of the solution. The comparison of truncation degrees 
n — 6 and n — 20 indicates that that the truncation degree n = 6 already yields a precise solution 
except when a2 is close to -1. The solution obtained by direct matrix inversion is preferable to 
the solution obtained by expansion in a2. Fig. ll2l shows the decrease in the values of the Legendre 
coefficients of the contraction solution up to ^ = 20 for various values of a2- The decrease is 
approximately exponential for \a2\ 1, but very slow for a2 — —1. The convergence of the 
harmonic expansion of the stress resultants is thus not guaranteed for extremal values of a2- 
Since the stresses and the strains include a supplementary factor a vanishing at the poles when 
a2 = —1, they are not affected by the convergence problem. 
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Figure 11: Influence of the truncation degree n on the Legendre coefficients of the contraction solution. The 
coefficients Ni (see Eq. ((T66J) are plotted as a function of 02 for £ = 0,2,4,6. The contraction is wo = —l and 
Poisson's ratio is equal to 0.25. The curve for £ = has been shifted, its real y-intercept being 2/(i/ — 1) = —8/3 
instead of —0.8. Three different approximations to the solution are plotted: n=6 with expansion limited to tSj 



(dotted curve). 



with no expansion in 02 (dashed curve), n=20 with no expansion in 02 (continuous curve). 
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Figure 12: Decrease with harmonic degree of the Legendre coefficients of the contraction solution (tDo = — 1). 
The coefficients Ni (see Eq. III66I 1) are plotted as a function of the harmonic degree (up to £ = 20) for different 
values of 02- The symbols overlap when £ = 0. Truncation degree n is equal to 50. 



Poisson's ratio has generally a weak effect on the solutions for contraction and despinning, 
except for the coefficients £ = for contraction and to a lesser extent i = 2 for despinning. The 
other coefficients are mainly affected by v at near extremal values of 0:2 ■ 

The computation of the stresses completes the example. The nondimensional stress resultants 
are computed with Eq. ([15]) and Eqs. d74t- ((76t : 

Nee = Fo + -F2 O2-P2 + -F4 + -Fe , (167) 
N^^ = FQ+F2 0iP2+FiOiPi + FGOiP6. (168) 

The Legendre coefficients of OiPj are given by the columns of the matrices ((80|) - (|8T|) . Ng^ is 
zero. The nondimensional stress is related to Nij by Eq. (IT5|) : 

a^J ^{l + a2P2) N,j. (169) 

The product of Legendre polynomials can be computed with the Adams-Neumann formula (jl43p 
but I will not give the explicit formulas. The stresses and the strains for the contraction solution 
at first order in 0:2 have already been given in Appendix 18.61 whereas those for the despinning 
solution are computed in Appendix 18.101 
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8.10 Despinning solution at first order 

If a depends only on the parameter a2 , the despinning stress function at first order in ck2 can be 
obtained from the contraction sohition (jl60p - (ll62p and the dualities (jl64p - (|165p : 

\ ( [2 21 + bv 6 I'-ll \\ , , 

P2 + a2\- + -— ^2 + :^-7^P4 . (170) 



W2 b + v\ \h 7 b + v 35 19 + j/ 

The despinning stresses are computed at first order in 0.2 from with Eqs. (IT5])- (fro|) : 

(5 + .)2(19 + .)- 







- 1 

— h 


W2 


5 + 




-< 




- 1 


W2 


54 





(172) 

'(5 + 1^)2(19 + :/) ^ ' 
where the auxiliary functions defined by 

see{'^,x) = - 31 + 18i/ + + (69 + 48j/ + 3i/2) a; + (40 + 18i^ + 2j/2) , (173) 

s^^(z/,x) = 21 - 30i/ - 31^2 + (47 + 72iy + t/2) X + (10 + 421/ + 8zy2) (174) 
have the property 

seeiv,l) = s^^{v,l) , (175) 

which ensures the equality of the stresses at the poles. In Eqs. (jl7ip - (|172p . P4 has been expressed 
in terms of P2'. 

A^|((P.)'-fp.4)^ <™' 

The strains can be computed at first order in 0.2 with Eqs. ([T7| - (fT9)) : 

^ {l-3iy)P2 + l + iy _ {I - ly^) eeeji^, P2) , . 

W2 b + u ' (5 + 1^)2(19 + 1/) ■ ^ ' 

^ {i~v)P2-l-u _ {l~y^)e^^{y,P2) 
W2 h + v ' (5 + z/)2(19 + z/) ■ ^ ' 

where the auxiliary functions defined by 

eee{v,x) ^ - 31 - 3i/ + (69 + i/) a; + (40 + 8i/) x^ , (179) 

e^^{v,x) = {2x + l){2l + v + (h + v)x) , (180) 

have the properties 

eg9{v,l) = eyy(i/, 1), (181) 

<egg{v,P2)> = 0, (182) 

e^^(z/,-l/2) = 0, (183) 

which ensure that the properties enumerated at the beginning of Section 13.11 are satisfied (note 
that < P2 1/4 and < (Fa)^ 11/32). 

At zeroth order in 0:2, the meridional strain is equal to its average value at = 45° (where 
P2 — 1/4), but this is not true at first order in 0:2 since eee{v, 1/4) ^0. If z/ = 0.25, the pseudo- 
nodes defined in Section [8^ occur at = 40.3° for the meridional strain (see Fig. [TJd), = 43.8° 
for the meridional stress and 9 — 64.4° for the azimuthal stress (see Fig. [SJd). 

Fig. 113b compares Eqs. (|17ip - (ll72p with the stresses computed at high order in 0:2, that is 
with truncation degree n = 20. The first order approximation is good for small 0.2 {0L2 — —2/5) 
but is rather bad for large 612 {a2 = —6/7). 
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(a) CONTRACTION (b) DESPINNING 
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Colatitude Idegreesl Colatitude IdegreesI 

Fi gure 13; Comparison of the stresses at first order in ck2 (dashed curves) and at truncation degree n = 20 
(continuous curves) : (a) contraction case (-100 = — 1), (b) despinning case («)2 = 1). The first order approximation 
is good if the equator-to-pole thickness ratio r = 1/2 (02 = ~2/5) and rather bad if r = 1/10 (02 = —6/7). 
Poisson's ratio is equal to 0.25. The meridional and azimuthal stresses are not distinguished but can be identified 
by looking at Fig. [5l 
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